

NASA Conference Publication 2199 

Spacecraft Dynamics 

as Related to 
Laboratory Experiments 

in Space 


Edited by 
G H. Fichtl 
B. N Antar 
and F G, Collins 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 


Proceedings of a workshop held at 
NASA Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 

May 1-2, 1979 


iwvsa 

Nat/onaf Aeronautics 
and Space Administration 


Sdantific and Teidinical 
Information Branch 





FOREWORD 


The Space Shuttle era will facilitate the conduct of 
scientific experiments in space on a routine basis. Many 
physics and chemistry laboratory experiments that will be 
performed in space will utilize the "low-gravity" feature 
of Earth orbit. However, the dynamic nature of spaceflight 
necessitates that scientists and engineers consider the 
effects of spacecraft dynamics in the design of experiments 
and planning of spaceflight missions. Prior to the meeting, 
the subject of spacecraft dynamics in the context of 
laboratory experiments in space hatd not been considered in 
a unified manner. The Physics and Chemistry Experiments 
in Space (PACE) Working Group sponsored this meeting for 
scientists and engineers to discuss the nature of spacecraft 
dynamics and its impact on the design and performance of 
laboratory experiments in space, and to identify engineering 
and scientific deficiencies in our knowledge of the subject. 
The meeting proved to be an effective forum. It is antici- 
pated that the proceedings will serve as an introduction to 
the subject for those new to the use of space as a laboratory 
environment and as a reference to those currently planning 
to use space for scientific experimentation. 

The support of Dr. Ellis Whiting, Program Manager for 
PACE, NASA Headquarters Office of Aeronautics and Space 
Technology, and the PACE Working Group for the successful 
completion of this workshop is gratefully acknowledged, as 
well as the participants for the preparation of their pre- 
sentations, manuscripts, and lively discussion. 
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INTRODUCTORY COMMENTS TO WORKSHOP ON SPACECRAFT 


DYNAMICS AS RELATED TO LABORATORY 
EXPERIMENTS IN SPACE 


During the 1980 's the Space Shuttle will provide a plat- 
form for scientists and technologists to perform experiments 
in space. In particular, the Space Shuttle will provide new 
opportunities for the performance of laboratory experiments 
in space which could not otherwise be performed on Earth 
because of the constraining effects of the Earth's gravi- 
tational field. However, the environment associated with 
Earth orbiting vehicles is not quiescent but, rather, is a 
dynamic environment which results from many sources, i.e., 
spacecraft maneuvers, vernier thrusters, atmospheric drag, 
crew activity, venting of liquids and gases, operation of 
machinery (blowers, fanS/ pumps, etc.) experiments, etc. 

(Table 1, p, 120) . The frequency bandwidth of this dynamic 
environment is broad, ranging from zero to frequencies asso- 
ciated with audible sound. It has been recognized that space- 
craft dynamics constitutes a major forcing function which 
could significantly affect the outcome of experiments in 
space and, as such, must be taken into account during the 
development, planning, operation, and postflight analysis of 
spaceflight experiments. However, the subject has not been 
fully discussed and explored in the context of spacecraft/ 
experiment interaction. This is particularly true for labo- 
ratory experiments to be performed in spacecraft that involve 
fluid dynamic processes, combustion processes, quantum fluids, 
crystal growth, life sciences, etc. The NASA Office of Aero- 
nautics and Space Technology (OAST) , Physics and Chemistry 
Experiments in Space (PACE) Working Group, has recognized the 
need for scientists, technologists, and efigineers to address 
the subject of spacecraft/experiment interaction. In this 
regard, the PACE Working Group has supported the conduct of 
this workshop relative to accomplishing the following objec- 
tives : 

1) To provide a forum wherein scientists and technolo- 
gists who are involved with the development and performance 
of experiments in space and specialists in the field of space- 
craft dynamics can exchange ideas with a view toward clari- 
fying the nature of spacecraft dynamics/experiment inter- 
action and developing understanding of the impact of space- 
craft dynamics on experiments in space. 
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2) To assess the knowledge that is available on the 
subject of spacecraft dynamics and to identify areas where 
important deficiencies exist in our information with respect 
to both spacecraft dynamics and the associated impact on 
experiments anticipated to be performed in space in the 1980 's. 

3) To develop recommendations on needed research and 
technology development to eliminate these deficiencies. 

Because of the broad nature of the subject of spacecraft 
dynamics/experiment interaction the meeting was restricted to 
laboratory experiments, e.g., free and contained fluids, 
quantum fluids, combustion processes, general relativity, and 
life sciences. The meeting was not concerned with Earth- or 
space-viewing experiments, except as related to viewing experi- 
ments' requirements which result in vehicle dynamics effects 
that might impact the laboratory experiments in question. The 
subsequent sections in this proceedings provide the papers 
presented at the meeting and deal with the following subjects: 

Part II - Low Gravity Experiments, Part III - Effects of 
Spacecraft Dynamics on Experiments, and Part IV - Spacecraft 
Dynamic Environment. A meeting agenda is provided in Appendix 
A. A list of participants in the workshop is available in 
Appendix B. 

During the first day of the workshop, a common frame of 
reference was established through presentations on laboratory 
experiments and spacecraft dynamics. Scientists and tech- 
nologists provided a brief description of laboratory experi- 
ments that are being planned for performance in space in the 
1980 's, including goals, objectives, and potential effects of 
spacecraft dynamics on these experiments. Presentations were 
also provided on effects of spacecraft dynamics on fluids in 
space by representatives from the scientific community. Subse- 
quent to these presentations, lectures were provided on space- 
craft dynamics encompassing orbital mechanics, rigid and 
elastic body dynamics, acoustic and vibration environments, 
and related subjects. In addition, on the morning of the 
second day two lectures were provided by Dr. Donald Drever of 
the California Institute of Technology on relativity experiments 
in space involving transmission of signals between satellites 
to detect gravitation waves, and by Mr, William Lange of the 
Bell Aerospace/Textron Corporation on spacecraft acceleration 
measurements with the Miniature Electrostatic Accelerometer 
(MESA) scheduled to fly on Spacelab Mission 2. The second day 
was devoted to open discussion aimed at accomplishing the pre- 
viously stated objectives. Scientific and engineering deficiencies 
were identified and the necessary developments needed to elimi- 
nate these deficiencies were identified. A summary of the 
deficiencies and needed developments is provided in the 
conclusions and recommendations section of Part I. 

George H. Fichtl 

Basil N. An tar 

Frank G. Collins 
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CONCLUSIONS AND RECOMMENDATIONS 


As a result of the presentations and discussions, impor- 
tant conclusions and recommendations were formulated. Major 
deficiencies were found to exist in our knowledge of space- 
craft dynamics and in the potential effect of spacecraft 
dynamics upon laboratory experiments in space. Recommenda- 
tions were made concerning technological and scientific 
developments that are needed to eliminate these deficiencies. 
These deficiencies and recommendations are summarized in this 
section. 

It was concluded by the workshop participants that there 
is an urgent need for a comprehensive description of the 
anticipated dynamic state of spacecraft for experiment defi- 
nition. This should include a complete description of the 
dynamic disturbance sources and strengths, and the descrip- 
tion should encompass a sufficiently wide bandwidth, i.e., 

0 to 50 Hz. For rotating fluid dynamic experiments wherein 
the angular velocity of a vessel is (rad sec“l) the 
description of the dynamics of spacecraft should encompass 
bandwidth 0 to Q/27 t (in units of Hertz) . Statistical 
descriptions of the dynamic environment of spacecraft 
(angular and rectilinear accelerations, angular rates, and 
attitude) are needed, including mean values, frequency 
spectra, discrete disturbance signatures as a function of 
risk of exceedance during the experiment duration time, dis- 
tribution functions, peak-to-peak values, etc. A major item 
upon which there was general agreement is the need for typi- 
cal examples of time histories and frequency spectra of 
spacecraft accelerations ( "g-jitter”) . It was further recom- 
mended that every effort be made by the NASA to acquire these 
profiles during early Shuttle missions (e.g., OFT missions, 
Spacelab 1 Mission) . Instrumentation should be provided on 
these early missions as well as all later missions to measure 
and record time histories of vehicle attitude, angular rates, 
and angular and rectilinear accelerations. Furthermore, 
event times of disturbances like thruster firings, gas 
venting, etc,, should be recorded for postflight analysis. 

In addition to the data above, information about high-fre- 
quency disturbances, e.g,, acoustic disturbances, is needed. 
The information described above is needed in general terms 
for experiment definition studies. Once a specific mission 
is identified for performance of the experiment, predictions 
of the dynamic environment for that mission are needed for 
experiment planning. Subsequent to flight, measurements of 
the dynamic environment are required for postflight experi- 
ment data analysis. 
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It was concluded that an assessment of the behavior of 
fluids in space in response to spacecraft dynamic behavior is 
required for experiment planning. The range of possible fluid 
f low/spececraft dynamics interactions that could occur is 
broad . Nevertheless , systematic development of a body of 
knowledge is possible, and it appears that it will be developed 
in part as a result of experiment payload definition and 
planning activities. It was pointed out that assessments of 
the role of spacecraft accelerations ("g-jitter”) on fluids 
may be required relative to destabilization of fluid experi- 
ments as a result of ”g- jitter,” effect of ”g-jitter" with 
zero mean value on convection, and masking of desired effects 
of ”g-jitter” deposition of energy into experiment systems. 

It was noted that in the case of experiments involving 
rotating fluids, important effects could occur from (1) pre- 
cession of the angular velocity vector of the experiment 
about the angular velocity vector of the orbit of the space- 
craft, (2) changes in vehicle attitude, and (3) periodicities 
in the rectilinear and angular accelerations of the vehicle 
at frequencies which are multiples of the magnitude of the 
angular velocity vector of the experiment. In the case of 
the latter, there appears to be essentially no information on 
the subject in the literature on rotating fluid mechanics. 
Finally, it was concluded that in certain cases the effect of 
spacecraft dynamics on experiments may not be predictable and 
could only be assessed by performing the spaceflight experi- 
ment. This is especially true of life sciences experiments. 

It was concluded that the following technological and 
scientific developments are required to help eliminate the 
deficiencies noted above. 

1) The available data on spacecraft dynamics from the 
Apollo, Skylab, and the Apollo-Soyuz missions should be ana- 
lyzed and documented. The results of these studies should be 
made available to the scientists who are planning experiments 
for space. 

2) Every attempt should be made to acquire data on 
spacecraft dynamics in the planned Shuttle flights in the 
1980 's. These data should include recordings from 3-axis 
low-g (10”9 g) accelerometers encompassing frequencies up to 
50 Hz. Also, to aid in experiment design it was recommended 
that VFI quality instruments to measure the dynatnic environ- 
ment should be flown on Shuttle flights. 

3) It was recommended that simulations be made of the 
dynamic environment of Spacelab using data from previous 
missions. Such simulations should take into account the 
forcing functions anticipated to be present during an experi- 
ment (e.g, , crew activity, venting, etc.), as well as our 


6 


i 



current knowledge of Orbiter systems (thruster time interval, 
specific impulse, limit cycles, etc.), 

4) It was concluded that improved understanding of the 
experiments can be made through more use of available ground- 
based facilities such as drop towers, laboratory simulation, 
computer simulation, intensive analytical models, and KC 135 
flights. It was concluded that information of this type 
would aid significantly in experiment design and planning, 

5) It was agreed that effort should be applied to the 
development of techniques for isolating experiments from 
spacecraft dynamic effects. 
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II. LOW-GRAVITY EXPERIMENTS 
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DEFINING CRITICAL POINT EXPERIMENTS 


FOR A SPACE LABORATORY 
by 

M. R. Moldover 
Thermophysics Division 
National Bureau of Standards 


ABSTRACT 


We are defining three representative low-gravity experi- 
ments for a fluid near its liquid-vapor critical point. Two 
of these experiments require very careful measurements of 
properties of the fluid in thermodynamic equilibrium, while 
the third experiment is a series of optical observations of 
the phenomena which occur as a fluid is changed from one 
phase to two phases, either by cooling through the critical 
point, or by adiabatic expansion. We are concerned with 
spacecraft dynamics insofar as residual spacecraft motions 
may complicate the interpretation of the data from the pro- 
posed experiments. It is possible that the spacelab environ- 
ment will render certain desirable experiments impractical. 


INTRODUCTION 


This experiment definition study has been sponsored by 
NASA Lewis Research Center under Interagency Agreement 
C-62861-C. An interim report has been prepared giving an 
overview of research in critical phenomena and phase tran- 
sitions in fluids with emphasis on the opportunities for 
experimental research that will become available in Space- 
lab [1] . A definitive review of gravity effects in fluids 
near the gas-liquid critical points has been published [2]. 

A brief overview and bibliography have also been published 
[3, 4]. 

A low-g environment greatly reduces certain constraints 
on the design of scientifically exciting experiments near 
critical points [1] . To establish the feasibility and 
utility of low-g experiments it is necessary to show that 
other likely constraints can be overcome. We are attempting 
to identify these constraints and methods of dealing with 
them for three representative experiments. 
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repeesentative equilibrium experiments 


Under normal gravity conditions, efforts to measure 
these equilibritom properties of fluids are ultimately frus- 
trated by the fluid's extraordinarily strong response to 
pressure gradients. This response is so strong near the 
critical point that the density at the bottom of a 1 mm high 
sample of fluid will be about 6% greater than the density at 
the top of the sample because of the sample's own weight. 

Thus any experiment which requires a sample at least 1 mm 
tall cannot yield meaningful results within ±3% of the 
critical density at the critical temperature. The density 
difference scales as g0»26. hence low-g experiments might 
yield meaningful results closer to the critical density. 
Residual steady spacecraft accelerations will prevent the 
experimental fluid samples from achieving the uniform 
density expected in the absence of such accelerations. 

The first equilibrium experiment we have studied in 
detail is a careful measurement of the angular dependence of 
the intensity and spectrum of light scattered from fluid 
near its critical point. This gives information about the 
size and lifetime of the transient "clusters" which occur 
near the critical point and are manifest as critical opales- 
cence. A schematic diagram of the experiment is shown in 
Figure 1. 

The second equilibrium experiment under consideration 
is a measurement of one of the theoretically important but 
experimentally "weak" anomalies which occur in several prop- 
erties of fluids near critical points. (The dielectric 
constant is our prototypical example of such a property; how- 
ever, the viscosity and constant volume specific heat, C^, 
are other properties of similar importance . ) A schematic 
diagram of a dielectric constant experiment is shown in 
Figure 2. A high-quality Earth-based dielectric constant 
measurement has been completed [5] . 

Under steady conditions near the critical point, a tem- 
perature gradient of 1 mK/cm will alter the fluid's density 
distribution about as much as an acceleration of 1 g. It is 
conceivable that unsteady spacecraft motion could deposit 
energy in fluid samples, leading to undesired temperature 
gradients. (Fluids near their critical points are extremely 
poor heat conductors.) A rate of energy deposition as small 
as 0.1 erg cm~3s"l would probably be a severe design con- 
straint for the proposed equilibrium measurements. Thus we 
need a method of estimating the rate of energy deposition in 
very compressible fluids samples from spacecraft motion. 
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Specific heat experiments, which require a complete 
record of energy inputs into a quite small fluid sample, 
will be particularly sensitive to unsteady motion. To 
establish the feasibility of a Gv experiment near the criti- 
cal point, consideration must be given to stray energy inputs 
to the fluid sample from mechanical vibrations of the system 
supporting the experiment and from acoustic sources in the 
air surrotinding the experiment. (Comparable Earth-based 
measurements are carried out in rather quiet laboratories.) 
Thus we require not only a thorough characterization of 
spacecraft dynamics, but also some sort of realistic upper 
bounds on the spectrum of noise and structural vibrations in 
Spacelab. 


NONEQUILIBRIUM EXPERIMENT 


A fluid near its critical density and temperature will 
be caused to phase separate by steady cooling and/or by 
increasing the volume of the cell several percent (see 
Figure 3). 

Under normal gravity conditions, the later stages of the 
phase separation experiment are dominated by sedimentation. 

Once bubbles (or droplets) are several microns in size, they rise 
(or fall) to the end of the cell within a few seconds. Under 
low-g conditions it should be possible to observe bubble and 
droplet growth for much longer periods of time. The growth 
regimes dominated by heat conduction and Marangoni convection 
should be more easily observed than on Earth. A detailed dis- 
cussion of the role of spacecraft dynamics on this phase 
separation experiment must await a more complete under- 
standing of the mechanisms dominating bubble motion under the 
nonequilibrium conditions of the proposed experiment. 
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DYNAMICS OF SUPERFLUID HELIUM 
IN ZERO GRAVITY 
by 

Peter V. Mason 

Low-Temperature Physics Group 
Jet Propulsion Laboratory 


Liquid helium becomes super fluid (Hell) when its temper- 
ature is reduced below 2.2 K. In this state, it has a number 
of unusual properties which are of technological and scien- 
tific interest. These include near-zero viscosity, thermal 
superconductivity, and the thermomechanical or fountain 
pressure. 

Hell is composed of two interpenetrating phases — the 
superfluid phase has zero viscosity, carries no heat, and has 
zero entropy. The other phase is a normal fluid with normal 
viscosity and entropy. The fraction of superfluid is deter- 
mined by temperature, rising rapidly from zero at 2.2 K to 1 
at 0 K. Hell has an extremely high thermal conductivity — 
about 1000 times that of copper — and can maintain temperature 
uniformity in a working volume of several centimeters 
uniform to a few milliKelvin. 

Hell may support several different sound waves, five or 
six at last count. First sound, the ordinary pressure- 
density wave of classical physics, is the synchronous motion 
of the two fluids. Second sound is the anti-synchronous 
motion of the two fluids without a change in density. It is 
a thermal wave; heat travels in a coherent manner described 
by a wave equation rather than the usual diffusion equation. 

Third soiind is still another form; in thin films the 
normal component is clamped to the substrate by its viscosity, 
but the superfluid, being frictionless, moves as a surface 
wave. In one gravity and at low frequency, this surface wave 
is a gravitational wave, similar to long- period ocean waves. 
High-frequency, short-period capillary waves, similar to wind- 
driven ripples on a water surface, are theoretically possible, 
but in actuality are so highly damped by energy exchange 
between the liquid and the gas above the surface that they 
cannot be observed. 

It is this capillary wave that we propose to observe in 
space. The low -gravity environment provides two essential 
experimental conditions. First, as shown in Figure 1, the 
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capillary wave behavior is extended to the long-period, low- 
frequency region, where attenuation is low enough to permit 
it to be observed. Secondly, the films must be about one 
micron thick over an area of several square centimeters, a 
condition impossible to obtain in the 1-g laboratory, due to 
tilts and irregularities in the substrate on which the films 
are formed. 

As shown in Figure 2, a time of flight technique will be 
used to measure the velocity and damping of the capillary 
waves. A pulsed heater will generate the wave, and a capaci- 
tance thickness detector will observe its passage. If the 
damping is low, the wave will pass over the sensors several 
times, permitting highly accurate determination of its 
velocity and amplitude loss. If the damping is high, only 
one passage will be observed, permitting determination but 
with reduced accuracy. 

Superfluid heliiam will be used as a cryogen to maintain 
the optics and detectors of space-borne far IR telescopes 
such as the free-flying Infrared Astronomical Satellite (due 
for launch in 1982) and the Shuttle Infrared Telescope 
Facility (planned for flight in the mid-to-late 1980 's). For 
these missions, it is important to understand the sloshing 
motion of the bulk liquid in near-zero gravity, for it may 
have deleterious effects on the sensors and the attitude 
control system. A preliminary 10-minute experiment has been 
performed on a rocket [1] supplemented by 15-second obser- 
vations in a NASA zero-g aircraft. A much longer and more 
accurate experiment will be flown on the Shuttle in con- 
junction with the surface wave experiment described above. 

As shown in Figure 3, the helium can be expected to pass 
from regions where surface tension dominates to regions where 
gravity dominates. Thus the behavior will be complex. In 
the surface tension region, the free surface will tend to 
take the shape of the vessel, while in the gravitational 
region, it will tend to be planar. Since the viscosity is 
low and nonlinear, theoretical predictions are difficult and 
direct observation is essential. The thruster and maneuver- 
ing systems of the Shuttle will provide the excitations of 
the liquid, and a large number of superconducting liquid 
vapor sensors will define the motion of the free surface. 

A third experiment will measure temperature distri- 
butions in the liquid and their correlations with the bulk 
motions. One expects variations of only millidegrees in 
volumes of the order of a cubic meter. We will use a number 
of sensors with a sensitivity of 10 microKelvin. 
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The overall system is shown in Figure 4. A helium cryo- 
stat will contain a sensor head with the experiments. A 
microcomputer will provide electrical excitations and gather 
and transmit the resultant responses to the experiment team 
on the groxind. Astronaut intervention will be possible, but 
will be used primarily for turn-on and shutdown, and for 
necessary commands when out of communcation with the experi- 
menters . 
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GEOPHYSICAL FLUID DYNAMICS 


by 

William W. Fowlis 

NASA/George C. Marshall Space Flight Center 


INTRODUCTION 


Geophysical fluid dynamics includes the study of large 
mass geophysical fluid flows, such as stellar circulations, 
planetary atmospheric circulations and ocean circulations. 
Systematic scaling or dimensional analysis reveals that cer- 
tain scales of these flows can be accurately modeled in the 
laboratory. Note that the procedure of laboratory geo- 
physical fluid flow modeling with which we are concerned is 
different from conventional engineering modeling. Rather 
than building a model to obtain numbers for a specific design 
problem, the relative effects of the significant forces are 
systematically varied in an attempt to deepen our under- 
standing of the effects of these forces. 

An area of geophysical fluid flow modeling that has 
received substantial attention is the modeling of large-scale 
planetary atmospheric flow in a rotating cylindrical annulus 
[1] . Figure 1 is a schematic diagram of the annulus. The 
boundary cylinders are maintained at different temperatures 
and the apparatus is mounted on a turntable. The volxxrae 
between the cylinders is filled with a low-viscosity liquid 
which has usually been water . Motions in the annulus have 
been made visible by adding aluminum powder to the liquid. 
Figure 2 is a photograph of an annulus flow taken with a 
camera rotating with the turntable and using a time exposure. 
An irregular jet flow showing a remarkable similarity to the 
atmospheric jet stream is observed. However, this similarity 
does not necessarily imply that the dynamics of the two flows 
are the same. Systematic scaling of the governing equations 
using characteristic values for both systems must be per- 
formed before we can conclude that the annulus is a model of 
the atmosphere. Such scaling of the horizontal momentum 
equation for the large-scale atmospherie motions and the 
annulus motions reveals a primary balance between the 
Coriolis force and the horizontal pressure gradient force for 
both systems. This is the geostrophic balance. A similar 
scaling of the vertical momentum equation reveals a hydro- 
static balance for both systems. A complete study of all the 
equations involved shows that the annulus is a model of the 
synoptic -scale atmospheric flow. 
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A major limitation of this annulus work is that it does 
not allow for the effects of spherical geometry. This is 
clearly important. For the annulus the rotation vector and 
gravity are always parallel, whereas for the Earth they vary 
from being parallel at the poles to being perpendicular at 
the equator. 


THE DIELECTRIC BODY FORCE 


There is no known way in which a true radial gravita- 
tional force can be obtained in a fixed geometry in the 
laboratory. However, a radial dielectric body force can be 
achieved. An analysis of the forces on a dielectric fluid in 
an electric field, E, reveals a dielectric force of the form, 
^E^Vg, due to variations in the dielectric constant [2] . For 
the geometry of a spherical capacitor (radii, Rj and Rq/ 

Rq > Ri) filled with a dielectric liquid and subjected to a 
voltage difference, V, across the spheres and thermal gra- 
dients on the boundaries, this dielectric body force, g^, is 
given by Equation (1) , where eg and pg are the dielectric 
constant and density 


^E 




i ^0”^I / 


( 1 ) 


of the liquid, respectively, at some reference temperature. 
The terms y and a are the thermal coefficients of dielectric 
constant and density, respectively; r is any radius between 
Rj and Rg [3] . 

A calculation for gg using realistic values reveals that 
it is not possible to make gg dominate terrestrial gravity. 
Thus, the need to perform these experiments in an orbiting 
laboratory like Spacelab becomes clear [4]. 


SPACELAB GEOPHYSICAL FLUID 
DYNAMICS EXPERIMENTS 


At the present time two spherical geophysical fluid 
dynamics experiments are being designed and prepared for 
Spacelab flights: 

1) The Geophysical Fluid Flow Cell (GFFC) , Principal 
Investigator: J. Hart, University of Colorado. 
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2) The Atmospheric General Circulation Experiment 
(AGCE) . Principal Investigator: W. Fowlis, NASA, Marshall 

Space Flight Center. 

It is anticipated that other geophysical fluid dynamics 
experiments will be proposed for later Spacelab flights. 

These two experiments are designed to examine different 
flows that have essentially different dynamics or physics. 

In the GFFC, the liquid between the spheres will be 
subjected to an unstable radial temperature gradient and 
rotation. The unstable gradient will be achieved by having 
the inner sphere everywhere wainner than the outer sphere. 

This is a model of stellar convection. For some experiments 
a latitudinal gradient will be imposed in addition to the 
unstable radial gradient. This may be a model of the Jovian 
atmospheric circulation. In the AGCE the liquid will be sub- 
jected to a stable radial gradient and an unstable latitudi- 
nal gradient and rotation. This is a model of the Earth's 
atmospheric circulation and extends the previous work per- 
formed in cylindrical geometry to spherical geometry. The 
thermal gradients will be imposed by maintaining temperature 
differences and gradients on the inner and outer spherical 
boundaries . 


MEASUREMENT AND DISTURBANCES 


A photochromic dye technique will be used for flow 
measurement [5] . When a small amount of a compatible photo- 
chromic dye CO. 01% by weight) has been added to the working 
fluid, it will darken upon exposure to ultraviolet radiation 
and then spontaneously clear again. Successive photographs 
showing the dye movement will reveal the flow. Figure 3 is a 
photograph showing photochromic dye sheets created in the 
working fluid in a spherical apparatus. This apparatus is 
being used for development work for the Spacelab experiments. 

A double-grid Schlieren system will be used for tempera- 
ture measurement [5] . In this system the image of a Ron chi 
ruling will be superposed back on the original ruling so that 
only light refracted by temperature gradients within the 
working fluid will pass through and be observed. For this 
system the inner sphere must behave like a mirror, and optics 
to correct for the spherical curvature are needed. 

Since these flows are thermally driven by temperature 
differences of about IOC® in confined volumes, flow rates 
will be small, 1 cm/sec and sometimes 0.1 cm/sec. Clearly, 



any disturbance of a rotational or precessional form to the 
rotating spherical containers could seriously affect the 
flows. 
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P - GENERAL POINT HAVING CYLINDRICAL POLAR 
COORDINATES (r. d, z) IN FRAME ROTATING 
WITH THE APPARATUS 

a = ( 0, 0, a) - ROTATION VECTOR 

g = CO, 0, -g) - ACCELERATION OF GRAVITY 

b, a, d - FLUID OCCUPIES REGION 

a< r < b, O < z < d [1 + a2(r2 - 1/2 (b2 + a2) ) / 2 gd] ^ d 

T (r, 0, z,t) - TEMPERATURE AT GENERAL POINT P AND TIME t 

Tg, Tb - T (a, 0, z. t); T ( b, 0, z. t ) RESPECTIVELY 


figure 1. SCHEMATIC DIAGRAM ILLUSTRATING ROTATING FLUID ANNULUS. 
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FIGURE 2. A TIME EXPOSURE SHOWING AN UPPER SURFACE WAVE FLOW 

IN THE ANNULUS. 
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FIGURE 3. PHOTOCHROMIC DYE SHEETS CREATED IN SPHERICAL APPARATUS 
BY EXPOSURE TO AN ULTRAVIOLET SOURCE. 
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DROP DYNAMICS 


by 

D . D . Elleman 
Jet Propulsion Lab 
Pasadena, California 


The Drop Dynamics Module is a Space lab-compatible 
acoustic positioning and control system for conducting drop 
dynamics experiments in space. It consists basically of a 
chamber, a drop injector system, an acoustic positioning 
system, and a data collection system. The principal means of 
collecting data is by a cinegraphic camera. The drop is 
positioned in the center of the chamber by forces created by 
standing acoustic waves generated in the nearly cubical 
chamber (about 12 cm on a side) . The drop can be spun or 
oscillated up to fission by varying the phase and amplitude 
of the acoustic waves. The system will be designed to per- 
form its experiments unattended, except for start-up and 
shutdown events and other unique events that require the 
attention of the Spacelab Payload Specialist. 

The proposed drop dynamics experiments will utilize the 
unique zero-g (low gravity) environment provided by space 
flight to investigate the dynamics of a free drop quanti- 
tatively for the first time. Aside from its fundamental 
academic interest, understanding of the behavior of a free 
drop will contribute to the sciences of nuclear physics, 
chemical processing, material processing and meteorology. 

This type of experiment is such a natural candidate for 
implementation in a manned Earth satellite that astronauts 
have already, on their own initiative, carried out limited 
qualitative experiments of this kind. The proposed experi- 
ments aim at obtaining precise data on the behavior of liquid 
drops by means of high-resolution cinematography in three 
orthogonal views. 

The theory of the dynamics of a free drop has been well 
studied in the approximation that dynamic quantities deviate 
linearly from a resting drop. With special exceptions to be 
discussed below, there is no nonlinear theory of the dynamics 
of a fluid drop. Not only are definitive experiments for the 
large amplitude behavior of fluid drops lacking but there is 
a deficit of definitive experiments even for linear behavior. 
This is a consequence of the limitations involved in con- 
ducting experiments in an Earth laboratory. Among these 
limitations are insufficient droplet sizes for accurate 
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observation, limited available time for experiments, and per- 
turbing effects due to the method of suspending the droplets. 

In the DDM a 1-cm radius drop will be positioned and 
held in an acoustic rotational well which can exert a force 
as large as 3 dynes on the drop. This is equivalent to the 
drops being able to withstand constant accelerations of the 
Spacelab of approximately 5 x 10“4 g. Acceleration forces 
larger than this will result in the drop striking the well of 
the chamber and a termination of that experiment. Dynamic 
accelerations smaller than the above level can also result in 
termination of the experiment and will be discussed. 
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INTERACTIONS BETWEEN SPACECRAFT MOTIONS AND THE 


ATMOSPHERIC CLOUD PHYSICS 
LABORATORY EXPERIMENTS 
by 

B. J. Anderson 

NASA/George C. Marshall Space Flight Center 


The Atmospheric Cloud Physics Laboratory (ACPL) is a set 
of equipment for studying a variety of atmospheric micro- 
physical processes in the low-gravity space environment. It 
will consist of three experiment chambers, an aerosol gener- 
ation and characterization system, and extensive support sub- 
systems. The experimentation planned for the first three 
ACPL missions will involve studies of the following processes 
the transport of heat, water vapor, and aerosol particles 
through air; cloud drop nucleation and growth; ice crystal 
growth from the vapor and within a cloud; drop freezing 
dynamics, aerosol production within clouds (SO 2 oxidation); 
and air turbulence decay. For each of these studies the 
advantage of working in space is a result of the absence of 
fluid motion in the experiment chamber, either fluid motion 
with respect to the chamber or with respect to the cloud 
particles themselves . On Earth these motions arise so easily 
from sedimentation and buoyancy forces that they are diffi- 
cult or impossible to suppress or eliminate. It follows then 
in the evaluation of the effects of spacecraft motions on 
ACPL experimentation that the motions of concern are those 
which will result in the movement of the fluid or cloud 
particles within the experiment chambers. 

Of the various vehicle motions and residual forces which 
can and will occur, three types appear most likely to damage 
the experimental results: non-steady rotations through a 

large angle, long-duration accelerations in a constant 
direction, and vibrations. Non-steady rotations present a 
problem because, for a fluid in a chamber, rotational motion 
of the chamber is poorly coupled to the fluid. The result is 
a long-duration rotation with respect to the chamber. The 
time constant for damping this motion is likely to be on the 
order of = 200 seconds for the ACPL expansion chamber. 

L is about 21 cm, j is the first Bessel function zero, and v 
is the kinematic viscosity. It appears that it is better to 
use the thrusters to stabilize the vehicle in a solid body 
rotation than to allow it to drift through large angular 
motions . 
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Constant accelerations, even of low magnitude, present a 
problem primarily because they may give rise to buoyancy 
forces when there are temperature gradients in the fluid. 
Apparently these forces do not fall off faster than /g. From 
a quick review of the literature it appears that the problem 
has not been adequately solved. If you calculate the Ray- 
leigh numbers for situations typical of ACPL (lengths of a 
few tens of centimeters, temperature differences on the 
order of 1®C, and 10"4g), you find the values are only 
slightly less than the critical Rayleigh number after cor- 
rection for a distributed heat source in the gas. Thus, it 
is possible to do the experiments without convection, but one 
must be careful. Some preliminary experiments may be called 
for to remove any doubt. Constant accelerations may also 
cause particle drift by sedimentation, but this concern is 
not critical for most ACPL work because cloud particles are 
generally quite small and viscous forces are significant. 

Concern over spacecraft vibrations arises in consider- 
ation of the ice crystal growth experinients. For most of the 
early ACPL work the crystals will be suspended near the end 
of a long fiber (20 cm long by 200 micron diameter) of glass 
or similar material. Clearly, small vibrations of the sup- 
ported end of the fiber could cause extensive motions of the 
ice crystal, if care is not taken to avoid this problem. 
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COMBUSTION EXPERIMENTS IN SPACE 


by 

A. L, Berlad 

State University of New York at Stony Brook 


Combustion phenomena of practical and fundamental inter- 
est occur primarily at normal gravitational conditions 
(g - 1) . The underlying natural convection-free (g = 0) com- 
bustion processes are obscured by gravitationally associated 
body forces. Accordingly, the Space Shuttle Laboratory could 
provide the laboratory conditions required for observation of 
fundamental combustion phenomena. These g = 0 observations, 
and associated analyses, are then expected to provide under- 
standing of combustion phenomena (at g ^ 0) which is not 
accessible by other means. A list of pivotal areas for 
needed combustion observations (available through Space 
Shuttle experimentation) includes: 

1) Single- and two-phase premixed flame propagation and 
extinction limits. 

2) Noncoherent flame propagation and extinction. 

3) Autoignition of premixed single-phase and two-phase 
combustible reactants . 

4) Upper pressure limit combustion phenomena and 
ignition, propagation and extinction processes in the 
neighborhood of such limits. 

5) Oscillatory combustion associated with the hydro- 
carbon-oxygen and with the carbon monoxide -oxygen systems. 

6) Two-phase flame spread and extinction phenomena 
involving large liquid-gas or solid-gas interfaces. 

7) Radiative ignition of solids and liquids. 

8) Pool burning. 

9) Smoldering of solid combustibles and the associated 
transition to flaming (or extinction) , 

10) Laminar gas jet combustion. 

11) Transient responses of combustible systems to time 
variations in gravitational field strengths. 
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A number of the preceding experiments may be impacted by unde- 
sirably high g- jitter effects. 

There are a number of compelling reasons, both experi- 
mental and theoretical, for taking combustion experiments 
into space. The first three that are listed in Figure 1 are 
general experimental eonsiderations and the last three are 
general theoretical considerations. When combustion experi- 
ments are conducted under normal gravitational conditions, 
large free convective forces tend to obscure the underlying 
combustion processes which are themselves determined by 
chemical kinetics, heat and mass and rtKjraentum transfer. The 
underlying molecular transport processes are thus generally 
clouded by these free convective processes. 

One of the really important experimental motivations for 
space-based experimentation is to assure the adequate/proper/ 
desirable preparation of the combustion experiment. For 
example, two-phase combustion experiments can be prepared in 
such a way that particulates don't settle. Accordingly, 
initial conditions for two-phase combustion experiments can 
be made uniform and homogeneous in space and time. Then, of 
course, the suppression of natural convection during the com- 
bustion process itself is made possible. Thus, one can 
characterize the space-based experimental observation of the 
combustion phenomenon as one which is not dominated by free 
convective forces, but rather one which is determined by the 
underlying chemical kinetics and molecular transport 
processes. From a theoretical point of view, this permits 
the resolution of a variety of key questions. These questions 
vary from experiment to experiment. Also, the theoretical 
representations of g = 1 combustion phenomena almost invari- 
ably become truncated or intractable when we try to include 
free convective processes at normal gravitational conditions. 
Space-based combustion experiments permit the use of real- 
istic theoretical representations which do not suffer from 
such problems. 

As a practical matter, most of us are interested in com- 
bustion behavior at normal gravitational conditions. Thus, a 
theoretical basis for g = 0 experimentation is that more 
meaningful, undistorted observation and more tractable theory 
can then be used to build on to understand the nonzero cases. 
This is the central feature of the space-based research in 
general, and, as indicated, detailed rationales vary from 
experiment to experiment. Twelve experiments are listed 
in Figure 2. Some of these we expect to be impacted by 
fluctuations in gravitational fields. 

Item (1) of Figure 2 refers to a class of experiments in 
which one looks at single- and two-phase flame propagation 
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and flame extinction. Here we do not expect to have signifi- 
cant effects due to some small perturbation in the gravi- 
tational field when g = 0. 

Item (2) refers to another class of experiments, so- 
called noncoherent flame propagation and extinction. This is 
an experiment for which one finds that rather than simply 
connected flame surfaces, the surface comes apart, flame 
propagation is noncoherent. Non-simply connected flame sur- 
faces are observed (at g = 1), and natural convection as well 
as molecular diffusion play a role in the experiment. Thus, 
the motivation for experiments of this kind would be to see 
whether in fact this phenomenon occurs at all and, if so, 
whether extinction limits are changed by the absence of a 
gravitational field. This phenomenon is particularly impor- 
tant at the so-called lean limit for H2~02. 

Autoignition or explosion is a criticality phenomenon 
(item 3) . One always worries about criticality phenomena if, 
in fact, we have perturbations in the neighborhood of a 
critical condition. 

The autoignition experiment for single- and two-phase 
systems is an experiment for which we have rather good g = 0 
theory. We do not have observations that are generally free 
of significant gravitational conditions. We expect that the 
observed autoignition condition would be affected by fluctu- 
ations in gravitational conditions. Heat loss mechanisms are 
very important in defining autoignition conditions, and where 
these are nonsteady or have average values which are somewhat 
different from what they ought to be at g == 0, this criti- 
cality condition can be affected. However, for g-values of 
the order of 10”4, substantial effects on this experiment are 
not expected. 

Item (4) is a very interesting experiment for which 
there is no current proposal. It is well known that at 
normal gravitational conditions and at very high pressures 
(of the order of 10-100 atmospheres) there is an upper 
pressure limit extinction phenomenon. It is difficult or 
impossible to propagate some flames under some conditions of 
high pressure. In fact, one of the effects of normal gravi- 
tational conditions is that natural convective and momentiam 
losses can be large, at high pressures. Extinction phenomena 
are generally associated with loss mechanisms. Thus, it 
would make a substantial difference on the operating body 
forces of these high-pressure, nonisothermal systems if we 
were to vary g- value over a substantial range, down to identi 
cally zero. There is no upper-pressure limit experiment 
currently proposed, though it would be a very interesting 
experiment indeed. 



Oscillatory combustion (item 5) is a very interesting 
and important theoretical and experimental issue in the field 
of combustion. We find oscillatory phenomena in the com- 
bustion of hydrocarbons with oxygen and in the oxidation of 
carbon monoxide with oxygen. The hydrocarbon-oxygen case is 
thought to be a largely thermokinetic oscillation, due to an 
interplay between a heat release rate and a heat transport 
rate through boundaries. The carbon monoxide -oxygen case is 
thought to be a purely kinetic phenomenon. This is perhaps 
the most isothermal of all nontrivial combustion phenomena. 

For the oscillatory combustion of carbon monoxide with oxygen, 
we find cycles of the order of several hundred. As many as 
450 have been observed in one experiment in which tempera- 
ture in the middle of this oscillating system was perhaps 
only 1°C warmer than that of the wall. Now, because mass 
transport of reactive species to walls is so important, one 
would like to experiment with some very large systems. How- 
ever, as we go to very large systems, we will run into natural 
convective processes. Although there appears to be an effect 
of physical scale and an effect of the chemical nature of 
walls g = 1, experiments of such a system are subject to free 
convective effects. 

Pool burning (item 6) is one combustion phenomenon where 
the effect of 1 to 10 Hz g- jitter may indeed be significant. 
Pool burning experimentation involves the propagation of a 
flame (that may be subject to extinction) over a stabilized 
pool of fuel. In such experiments, the ability to establish 
a very well defined initial condition and maintain it during 
the course of the experiment is important. 

Item (7) is concerned with two-phase flame spread 
and extinction phenomena. Such an experiment may involve a 
matrix-like structure of particulates, with flame spreading 
occurring over the particulates. Such an experiment has been 
identified previously as one of significant interest from 
both fundamental and applied points of view. 

Item (8) , radiative ignition of solids and liquids, is 
thought to be an important experiment, to be done at reduced 
gravitation or zero gravity fields. 

Smoldering (item 9) is perhaps one of the most important 
phenomena from the point of view of fire safety at g = 1. 
Combustion characteristics of cigarette, cigar, and mattress 
burning processes involve transitions of smoldering to 
extinction, smoldering to flaming, or just continued 
smoldering. Smoldering itself is not well understood theo- 
retically. Neither are smoldering-related transitions to 
flaming or extinction. An experiment which is designed to 
study the transition of smoldering to either flaming or 
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extinction which is going to develop over a period of time 
could, in fact, be significantly influenced by whether or not 
we have a gravitational field. All of the processes which 
participate in the smoldering phenomenon appear to be rather 
slow. Mass and heat transport as well as chemical kinetic 
and heat release rates are all slow. Natural convective 
processes, though modest, may still be very important for 
this class of phenomena. 

Laminar gas jet combustion (item 10) studies have 
already been conducted in drop tower experiments. Many of 
these studies involve rather substantial forced convective 
conditions and we don't expect that small gravitational per- 
turbations will create significant effects. 

Criticality conditions of various kinds are common in 
combustion phenomena. There are flames that do exist; others 
don't exist; there are transitions to extinctions; there are 
transitions to detonation from ordinary flame propagation; 
and one would expect that those transition conditions may be 
affected if time-dependent gravitational conditions are 
present. This is indicated for the class of item (12). 
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COMBUSTION EXPERIMENTS IN SPACE 

OBSERVATION AT g O PERMITS 

(1) PROPER/DESIRED PREPARATION OF EXPERIMENT. 

(2) SUPPRESSION OF NATURAL CONVECTION - PRIOR TO & DURING COMBUSTION. 

(3) CHARACTERIZATION OF UNDERLYING COMBUSTION PHENOMENON. 

(4) RESOLUTION OF A CENTRAL QUESTION IN COMBUSTION THEORY/EXPERIMENT. 

(5) MORE REALISTIC - MORE TRACTABLE THEORY. 

(6) THEORETICAL BASES FOR g > O. 


FIGURE (1) 
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EXPERIMENTS 


(1) SINGLE/TWO-PHASE FLAME PROPAGATION AND EXTINCTION. 

(2) NONCOHERENT FLAME PROPAGATION AND EXTINCTION. 

*{3) AUTOIGNITION - SINGLE/TWO-PHASE. 

*(4) UPPER PRESSURE LIMITS - IGNITION, EXTINCTION, PROPAGATION. 

*(5) OSCILLATORY COMBUSTION: 

[HYDROCARBON - OXYGEN] + [CARBON MONOXIDE - OXYGEN] 
*(6) POOL BURNING. 

*(7) TWO-PHASE FLAME SPREAD AND EXTINCTION. 

(8) RADIATIVE IGNITION OF SOLIDS AND LIQUIDS. 

*(9) SMOLDERING OF SOLID COMBUSTIBLES AND TRANSITIONS TO 

FLAMING/EXTINCTION. 

(10) LAMINAR GAS JET COMBUSTION. 

*(11) 0<g«1. LOW (g) STUDIES. 

*(12) TRANSIENT RESPONSES OF COMBUSTIBLE SYSTEMS TO TIME VARIATIONS IN 
GRAVITATIONAL FIELD STRENGTHS. 

*POSSIBLE g-GITTER EFFECTS OF SIGNIFICANCE. 


FIGURE (2) 
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TWOr-PHASE FLOW AND HEAT TRANSFER 


UNDER LOW GRAVITY 
by 

Walter Frost 

The University of Tennessee Space Institute 
Tul lahoma , Tennes see 


The two-phase flow and heat transfer program for the 
space laboratory in which I am involved is outlined in 
Figure 1. General Dynamics is under eontract with NASA, 

Lewis Research Center, to design a two-phase fluid mechanics 
heat transfer facility for spacecraft. I *m listed as the 
lead scientist for the two-phase flow patterns, pressure 
drop, and flow boiling experiment. Three other projects 
being considered under this same conceptual design study are 
listed in Figure 2; liquid reorientation, pool boiling, and 
bubble dynamics. Only liquid-vapor, two-phase flows are 
considered in this paper. 

The needs for the spacecraft experiment described above 
are twofold. One is for space transportation systems which 
will begin to develop in the 1980 's. A number of these 
systems will involve processes containing two-phase fluid 
flow and heat transfer. General Dynamics in their previous 
study II] has identified the five areas shown on Figure 3 
for which some 45 different applications where two-phase flow 
will occur have been identified. An example of one common 
occurrence of two-phase flow and heat transfer in these 
applications is cryogenic transfer lines where a cryogen is 
pumped into a warm line and undergoes boiling until the line 
cools to the appropriate temperature. Power generation 
systems utilizing boilers and condensers in which two-phase 
flow and heat transfer occur are also an example of typical 
applications . Since two-phase phenomena occur in these 
applications for a variety of fluids ranging from hydrogen to 
liquid metals, two-phase flow and heat transfer predictive 
models are needed for the design of space transportation 
systems . 

Current design correlations for predicting two-phase 
flow pressure drops and heat transfer are primarily empirical 
and have been in most all cases developed under standard 
gravity. The flow process is generally so complex that 
writing the exact governing equations is prohibited and conse- 
quently we can’t look fundamentally at how variation in the 
dependent variable will affect the heat transfer process. 
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Thus, there is a need to carry out experiments under reduced 
gravity to determine how it influences these correlations. 

A second need to study two-phase flow and heat transfer 
in space is to improve our fundamental understanding of this 
phenomenon. By isolating gravity forces from surface tension 
forces, inertia forces, and other forces that are involved in 
the flow and boiling process, we will be able to establish a 
meaningful theoretical model to replace the empirical tech- 
niques currently employed. 

Thus, there is a twofold need for two-phase fluid 
experiments in space. One is to develop design criteria that 
can be used for future space system design, without neces- 
sarily understanding completely the physical mechanism 
taking place, i.e., empirical design correlations. The other 
is to carry out fundamental experiments to provide a better 
understanding of the physics of the processes for future 
analytical work. 

Although there are several classifications in boiling 
[2], we'll skip these and begin by considering in detail the 
two-phase flow pattern which occurs in a heated vertical 
duct. Figure 4. Fluid is pumped upward through the duct. 

The duct is uniformly heated along its length, and as fluid is 
pumped upward it changes phases from pure liquid to pure 
vapor. 

On the left of the figure the temperature variation 
along the tube length is plotted. Both the wall temperature 
of the duct and the bulk temperature of the fluid in the duct 
are plotted. The visual flow patterns which occur are listed 
on the right-hand side of the figure. The fluid enters as a 
pure liquid, and heating takes place initially by standard 
convective heat transfer. Thus, single-phase flow initially 
takes place in the entrance region , Water is considered in 
our discussion, but all fluids will behave similarly. As the 
thermal boundary layer builds up on the duct wall and the 
temperature exceeds the saturation temperature of the fluid, 
bubbles begin to form at specif ie nuclei on the surface. 

These bubbles will grow, but once they grow outside the 
thermal layer, they collapse because the bulk fluid tempera- 
ture outside the thermal layer is below the saturation 
temperature of the fluid — the dashed line in the figure. 

This form of boiling is called subcooled nucleate boiling 
because the bulk fluid temperature is below the saturation 
temperature. It is also referred to as surface or local 
boiling because the bubbles exist only in the location of 
the surface and not within the fluid stream; i.e,, they 
collapse once outside the thermal layer required to make them 
grow. However, as the flow continues upwards and energy is 



continuously added, the bulk temperature of the liquid 
reaches the saturation temperature of the fluid. Now the 
bubbles will continue to grow even outside the thermal 
boundary layer. The flow takes on a bubbly or frothy flow 
pattern . 

As the fluid moves upward and energy is continuously 
added, nucleate boiling continues to take place on the heated 
duct surface; however, the bubbles become sufficiently large 
that they now coalesce and form large slugs or a frothy 
mixture of vapor. This pattern is called the slug, churn, 
or froth flow pattern. 

Soon the vapor generation becomes so large, and the 
vapor flow so fast that it forces the liquid to the walls of 
the duct and the vapor flows up the center of the duct as a 
continuous core. This is called the annular flow pattern. 
Further downstream an exchange of energy between the vapor 
and the liquid layer takes place through droplets being torn 
from the surface of the liquid and then reentrained into the 
layer downstream. This is called the spray or liquid dis- 
perse pattern, or sometimes the annular mist region because 
the liquid annulus exists with mist in the vapor core. 

Eventually the liquid layer will dry out, which is 
called the point of dry-out or burn-out. The wall tempera- 
ture jumps drastically when the liquid layer dries out 
because a very effective mechanism of cooling exists as long 
as liquid covers the surface. Once this breaks down, the 
heat transfer coefficient to the pure vapor is comparatively 
low and the v/all temperature rises rapidly. The flow pattern 
then becomes a mist with droplets in the entrained vapor, and, 
finally, all liquid evaporates and single-phase vapor flow 
occurs. 

All of the previously described flow patterns are visual 
definitions of the flow process. The difficulty with pre- 
dicting flow and heat transfer in these various flow pattern 
regimes is that associated with the different flow patterns 
are different mechanisms of pressure drops and heat transfer. 
Convective heat transfer to a pure liquid or vapor is the 
mechanism which occurs in the single-phase regime; however, 
we will not discuss this mechanism. Rather, consider the 
nucleate boiling process which occurs on the walls in the 
subcooled boiling and slug flow regimes. Nucleate boiling 
heat transfer is a very efficient heat transfer mechanism. 

The bubbles growing, collapsing, and being detached from the 
surface create a very high transport of energy next to the 
wall. Arguments still exist, however, as to the exact 
mechanism by which the bubbles enhance the heat transfer. 

The contribution due to latent heat of evaporation as 
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compared to that due to augmented turbulence effects is not 
clearly defined. Although a few other factors are involved, 
we know that the bubbles growing and collapsing are 
responsible for the very high heat transfer rates obtained. 

Nucleate boiling continues on the wall through the slug 
flow regime. The mechanism of nucleate boiling heat transfer 
has been well studied, and the influence of a number of 
variables dociamented 12,3] . Equations governing the bubble 
dynamics have been written, and how gravity enters into the 
particular growth rate and transient forces acting on the 
bubbles can be estimated. Near the wall the bubble growth is 
basically due to thermal effects; i.e., it depends upon the 
rate at which energy is transported to the liquid-vapor 
interface to evaporate the liquid. The bubble growth is, 
thus, pretty much independent of any gravity forces during 
this stage of growth. This has been partially confirmed by drop 
tower and flight tests, at least for boiling in a large pool 
of liquid. Natural convection occurs in large pools rather 
than forced convection, as illustrated in Figure 4 . 

The characteristics of pool boiling heat transfer have 
been outlined by Siegel [3] as shown in Figure 5. In the 
figure heat flux is plotted versus temperature difference 
between the wall and the saturation temperature of the fluid. 

The curve is essentially independent of gravity, although 
very small gravity effects have been observed in drop tower 
tests or flight tests. These are generally negligible for 
design purposes. Thus, gravity is not expected to have a 
pronounced effect on the mechanism of nucleate boiling heat 
transfer near the wall. 

Gravity will, however, influence the size of bubbles 
which are torn off the wall by other forces. The bubble 
departure from the wall is dependent upon four basic forces; 
one is the surface tension force which tends to hold the 
bubble on the wall; another is the viscous forces which in 
the forced flow case tend to drag the bubble off the wall— in 
a pool boiling experiment they tend to hold the bubble on the 
wall; the buoyancy or gravity force which lifts the bubble 
off the wall; and the inertia force which is a function of 
how fast the bubble grows— if it grows very rapidly, it causes 
the liquid to have a certain amount of inertia, and this 
inertia can force the bubble away from the wall. The 
competing forces of gravity and of inertia are basically the 
forces influencing the size at which the bubble is torn off 
the wall. Certain types of bubbles are strongly dependent 
upon gravity and for low gravity these bubbles will become 
very large before they are lifted off the wall. Consequently, 
slugs and churn flow will probably occur quicker at low 
gravity. 



The effect of gravity on isolated bubbles in a fluid has 
also been studied. We know obviously that bubbles will rise 
faster under high gravity. Once again, however, there will be 
competing forces which are not fully understood. For example, 
the drag versus gravity effects. Utilizing experiments in 
space, the influence of drag forces can be isolated from 
gravity. 

In the nucleate boiling regime, we expect that the sub- 
cooled transfer mechanism is not strongly influenced by 
gravity. However, the transition from nucleate boiling to 
slug flow will probably occur quicker because it is dependent 
upon the rate the bubbles rise and coalesce with one another 
and form slugs. 

Since the bubbles cannot escape as fast at low gravity, 
the occurrence of slug flow is enhanced. There are also 
competing forces in the transition from slug to annular flow. 
This transition is dependent on a force balance between 
liquid flowing down the wall by gravity and vapor dragging 
the fluid up the wall. Hence, the important dry-out or 
burn-out conditions will be strongly influenced by gravity. 
Thus, there is a lot to be learned about the physics of two- 
phase flow and heat transfer by carrying out an experiment in 
space. 

In the annular flow region the mechanism of heat trans- 
fer is not nucleate boiling. Rather it occurs by conduction 
or convection through the thin liquid layer and evaporation 
at the interface. Again, this is a very effective mechanism 
for transferring heat, but all present design correlations 
are based on experimental results under earth’s gravity, 

A second complicated problem is that although we can 
approximately calculate pressure drop and heat transfer in 
the various flow regions described above, we need to know 
under what conditions these regions occur in order to employ 
the appropriate design criteria. To date, the conditions 
under which these regions exist are not very well defined. 
Experiments have been carried out using air and water 
mixtures. In general the liquid flow rate is held constant 
and vapor flow rate is increased systematically. Figure 6 
illustrates a flow regime map measured by Griffith [4] , 
whereas the dashed lines are regime boundaries recommended 
by Kazlov 15] . The flow regimes are bubbly flow, slug and 
froth flow, and mist annular flow as simulated using air- 
liquid mixtures. The gravity dependence of the flow regimes 
is based on dimensional analysis and has not been confirmed 
with experiments at sustained low gravity. 



This influence of gravity may be even more significant 
for horizontal flows because now the fluid can stratify, as 
illustrated in Figure 7 , Increasing vapor flow with the 
liquid flow constant results in the different flow patterns 
illustrated in the figure. These particular flow conditions 
may vanish altogether in a reduced gravity environment and, 
in fact, there is evidence supporting this from flight 
experiments conducted by General Dynamics II] , 

Current techniques available to predict the flow regime 
are thus in the form of flow regime maps , Figure 6 , Flow 
regime maps are based on earth gravity controlled experiments. 

It should be noted that the lines on the maps are really 
broad bands representing transitions between the different 
flow regimes. The empirical results obtained under earth's 
gravity where g is introduced by dimensional reasoning indi- 
cate that as gravity becomes very low, slug and froth flow 
may not occur, and that annular flow will be the dominant 
flow regime. Thus, we expect that slug and froth flow might 
not exist in space. However, we need to verify this. 

The purpose of the Spacelab experiment is to do three 
things directed toward this goal. First, investigate two- 
phase flow patterns under reduced gravity using the water-air 
mixture experiment. A schematic of the proposed test 
apparatus is shown in Figure 8. Air and water are circulated 
through the system. The quality or the mixture of air-water 
is controlled. Photographs of the test section will be made 
and at the same time pressure drops across the test section 
will be measured. The data will establish a flow regime map 
under reduced gravity conditions with corresponding pressure 
drop correlations. 

The test section is also equipped with an electrical 
resistance heater in order to allow a flow boiling experiment 
to be carried out; Freon 11 will be used. Again, high-speed 
photographs will be taken of the test section to determine 
flow patterns. Measurements of the temperature gradient and 
pressure drop along the duct will also be made. Thus, 
quality change can be measured, and heat transfer calculated. 
This will provide information on two-phase flow and heat transfer 
under low-gravity conditions. More details of the experi- 
mental program are given in Reference 6, 
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OTHER PROJECTS AND LEAD SCIENTI ST 


• LIQUID REORIENTATION 
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FIGURE 2 OTHER PROJECTS AND LEAD SCIENTIST. 
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• DESIGN AND PERFORMANCE OF SPACE TRANSPORTATION SYSTEMS 


• PROPULSION 

• AUTOMATED SPACECRAFT 

• POWER GENERATION 

• EXPERIMENTAL SUPPORT (SPACELAB AND SPACE STATION) 

• ENVIRONMENTAL CONTROL/LIFE SUPPORT SYSTEMS (EC/LSS) 

• ADVANCE FUNDAMENTAL UNDERSTANDING OF TWO-PHASE FLOW 
AND BOILING HEAT TRANSFER 


FIGURE 3 NEED FOR EXPERIMENT 
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FIGURE 4 CHARACTERISTICS OF FORCED CONVECTION TWO-PHASE HEAT TRANSFER. 
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1. FOR NUCLEATE BOILING, BOTH ANALYSIS AND EXPERIMENT INDICATE THAT 
THE RELATION BETWEEN TEMPERATURE DIFFERENCE (Tw - TSAT ) AND 
WALL HEAT FLUX IS INSENSITIVE TO GRAVITY. 

2. THE PEAK NUCLEATE BOILING (CRITICAL) HEAT FLUX WAS FOUND EXPERI- 
MENTALLY TO VARY REASONABLY WELL AS G’A, WHICH IS IN AGREEMENT 
WITH THEORY. 

3. THE RELATION BETWEEN Tyv - TsAT AND WALL HEAT FLUX IN THE TRANSI- 
TION REGION BETWEEN NUCLEATE AND FILM BOILING APPEARS FROM 
LIMITED DATA TO BE INSENSITIVE TO GRAVITY REDUCTIONS. 

4. THE MINIMUM HEAT FLUX VALUE WHERE TRANSITION BOILING CHANGES TO 
FILM BOILING DEPENDS ON G%. 

5. IN LAMINAR FILM BOI LING THE HEAT-TRANSFER COEFFICIENT DEPENDS ON 
G%. FOR A TURBULENT VAPOR FILM THE EXPONENT INCREASES AND MAY BE 
AS LARGE AS 2/5 TO 1/2. 


FIGURE 5 CHARACTERISTIC OF POOL BOILING HEAT TRANSFER [31, 
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FIGURE 7 FLOW PATTERNS IN HORIZONTAL TUBES 
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SPACELAB MODULE 



FIGURE 8a FLOW PAnERN AND PRESSURE DROP EXPER|jV\ENT FLOW SCHEMATIC 
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TRIBOLOGY EXPERIMENT 


by 

William A. Wall 

NASA/ George C. Marshall Space Flight Center 


The purpose of the Tribological Experiments in Zero 
Gravity (TEZG) , Figure 1, is to examine the interaction of 
liquid lubricants and surfaces under static and dynamic con- 
ditions in a low-gravity environment. Three functional 
objectives are planned; (1) Fluid Wetting and Spreading 
Experiments (FWS) ; (2) Journal Bearing Experiments (JB) ; and 

(3) means to accurately measure and record the low-gravity 
environment during experimentation. The scientific objective 
of the fluid wetting and spreading portion of the TEZG experi- 
ment is to study fluid-surface interactions lander static con- 
ditions. This is to be done by observing the wetting and 
spreading process of selected commercial lubricants on repre- 
sentative surfaces in a near- zero gravity environment. 

The useful life of a machine element is often determined 
by the slow migration of a liquid lubricant on its surface. 
This phenomenon is commonly known as wetting and spreading 
and is controlled by the minute interfacial forces of the 
three material phases; namely, the solid surface of the 
machine element, the layer of the liquid lubricant, and the 
surrounding air. When a static state prevails, the boundary 
of the lubricant coverage is stationary; then the situation 
can be characterized by the contact angle together with the 
surface tension of the lubricant. However, during a wetting 
and spreading process, the contact angle theoretically reduces 
to zero, and the rate at which the front of the lubricant 
advances can no longer be completely described by the static 
properties. Accurate observation of the wetting and 
spreading process in an earth environment is not possible 
because earth gravity is many orders of magnitude larger than 
the surface forces and would thus dominate the fluid motion. 
The prolonged zero-g condition during a Spacelab mission pro- 
vides a unique opportunity to acquire quantitative knowledge 
of the wetting and spreading phenomenon. 

Since FWS concerns both a fluid and a surface, it is of 
interest to include in the experiment several fluid-surface 
combinations. The present plan calls for grouping numerous 
surface coupons into an assembly. Various fluids will be 
brought to the surface through a small centrally located duct. 
An automated dispensing mechanism permits the process to 
start simultaneously on all coupons of the assembly. Also, a 
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photographic record will readily depict any significant 
differences in the wetting and spreading rates. 

The scientific objective of the journal bearing (JB) 
portion of the TEZG experiment is to study hydrodynamic films 
in journal bearings operating in zero gravity and the 
stability of the dynamic bearing system under these con- 
ditions. Efforts will be made to control the development of 
the fluid film by use of experimental bearing configurations. 
Under dynamic conditions, such as in a simple journal bearing, 
the formation of hydrodynamic films is affected by low 
gravity. Because the ability of a journal bearing to develop 
a hydrodynamic, load-supporting film is responsible for the 
universal use of such bearing systems, these fluid surface 
reactions are very important. The film thickness developed 
in such a bearing is a fxinction of the surface speed N, the 
Ixabricant viscosity u, and an inverse function of the load W 


h = f 


Nu 

W 


In near- zero gravity, the load W approaches zero and the 
bearing system should go unstable. Data on the development 
of instability xinder low-gravity conditions will be of great 
interest to scientists in the fields of bearings and Ixibri- 
cation. Since, in theory, a nonloaded rotating journal is not 
stable, it is planned to investigate the possibility of using 
experimental bearing configurations to control the development 
and location of the hydrodynamic film. Results of these 
tests should be of great interest to hydrodynamic is ts as the 
development of these data is impossible except under low- 
gravity conditions. It is believed that the information 
obtained during these experiments will provide a wealth of 
data on the development of hydrodynamic films and journal 
bearing stability under no-load conditions. 

To verify the accelerations under which these experi- 
ments are subjected, a 3-axis, low-gravity accelerometer will 
be located in the TEZG experiments operations drawer adjacent 
to the experiments during operation. This accelerometer 
system has a threshold of 1 x 10”4 g»s and a maximum output 
of 31 X 10~3 g's. It is desired that the TEZG experiments be 
run in an environment of less than 1 x 10“3 g's for a success- 
ful experiment. 
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Figure 1. Tribology Experiment Assembly 



III. EFFECTS OF SPACECRAFT DYNAMICS ON EXPERIMENTS 
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SUSCEPTIBILITY OP MATERIALS PROCESSING EXPERIMENTS 
TO LOW-LEVEL ACCELERATIONS 


by 

Robert J . Naumann 

NASA/George C. Marshall Space Flight Center 


Recognizing that it is impractical to maintain an actual 
zero-gravity environment in an orbiting vehicle,! the 
question is frequently asked, how much residual acceleration 
can an experiment tolerate? This is not a simple question to 
answer for the reasons that, (1) there are a multitude of 
different types of low-level disturbances, and (2) there are 
a multitude of types of experiments that have varying degrees 
of susceptibility to each type of disturbance. 

In order to attempt to answer this question, it is neces- 
sary first to enumerate the types of disturbances expected. 
These may be classified into four categories; (1) steady 
low-level accelerations, (2) compensated transient acceler- 
ations, (3) uncompensated transient accelerations, and 
(4) rotations. The origins and effects of the disturbances 
will be examined briefly in the following paragraphs. 

By steady low-level accelerations, we mean accelerations 
on the order of 10"5 g and smaller that are relatively con- 
stant for durations of tens of minutes or longer. These 
arise from aerodynamic drag, Keplerian effects, and slow 
vehicle rotations. Accelerations arising from aerodynamic 
drag, of course, depend on vehicle altitude, orientation, and 
atmospheric density which has a day-night variation and 
depends to a lesser degree on solar activity. This force is 
constant in direction for a vehicle in earth-fixed orien- 
tation. For a vehicle in an inertial orientation this force 
rotates in the orbital plane with the orbital period, 
resulting in a nearly zero time average over one orbital 
period . 

Keplerian forces originate because every particle in the 
spacecraft has its own individual orbit. These orbits in 
general have slightly different periods and would result in 
oscillations of the particles coupled with a drift along the 
orbital path, unless they are constrained. The magnitude of 


1. This could only be accomplished by the heroic measure of 
actually flying the vehicle around the experiment in question. 
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the constraining force required is on the order of 10 '-'7 g per 
meter of displacement from the center of mass of the vehicle. 
The direction of these forces generally oscillates with a 
period of twice the orbital period. 

Low-level, nearly steady accelerations are also produced 
by slow vehicle rotations such as those produced by the 
rotation required to maintain earth-fixed orientation, those 
produced by gravity gradient torques in a freer^drift mode 
with altitude control thrusters shut off, or those produced 
by the barbecue roll sometimes required for thermal control , 

The magnitude of this acceleration ranges from 10'~7 g per 
meter displacement from the center of mass in the case of one 
revolution per orbit, to approximately 10”6 g per meter for 
the other types of rotation. In all cases the direction of 
this force is along the radius vector from the vehicle center 
of mass to the experiment in question. 

Grashof numbers can be as much as ~20 for a steady 
acceleration of 10~6 g^ assuming a container with character- 
istic dimensions of 10 cm containing H2O with a 10° tempera- 
ture differential. This would result in a flow of ~40 lam/sec. 
Such flows are sufficient to traverse the container in one- 
half the orbital period. Therefore, considerable mixing can 
occur in an experiment even if conditions are chosen to 
average out low-level accelerations over an orbital period. 

By compensated transient accelerations we mean short- 
duration impulses that dd not involve a momentum change to 
the spacecraft. Such transient aceelerations arise from crew 
activity or internal mechanical vibrations. The important 
distinction is that such impulses are compensated a short 
time later by equal and opposite impulses. Even though such 
impulses may be as high as 10~2 g^ the flow velocity acts for 
only a short time and the displacement is small. Furthermore, 
the random nature of these disturbances further reduces the net 
displacement. Therefore the fluid transport resulting from g- 
jitter (as such motion is often referred) is negligible for 
most cases. However, such transients may be of concern in 
containerless suspension systems or in experiments with free 
surfaces such as liquid floating zones. 

Uncompensated transient accelerations result from 
external forces applied to the spacecraft such as thruster 
firings for attitude control or orbital maneuvers. In Skylab, 
attitude control was generally accomplished by means of 
control moment gyros which did not impart external forces to 
the spacecraft. It was necessary, however, to periodically 
dump the angular momentum accumulated from the integrated 
gravity gradient torques by means of cold gas thrusters. 



Also, roll maneuvers to execute various earth resource tasks 
were accomplished by means of thrusters . In the Apollo-Soyuz 
mission, attitude control thrusters fired as a couple which 
exerts a torque on the vehicle without a change in linear 
momentvim. Therefore, no uncompensated transient accelera- 
tions were experienced. 

Shuttle, on the other hand, uses only downward firing 
thrusters for pitch and roll control. This will result in 
uncompensated transient accelerations with each attitude 
control pulse which has a magnitude of 10'"4 g. These acceler- 
ations can produce fluid velocities as high as 4 mm/sec in 
the 10 cm cell described previously. However, since the 
duration of these pulses is only 100 milliseconds, the 
resultant displacement is on the order of 100 microns for 
each thruster firing. Since these thrusters always fire in 
the same direction relative to the vehicle, the resultant 
displacement will be cumulative. The net result for free 
convection is similar to the case of a steady acceleration 
with a magnitude equal to the time averaged value of the 
thruster accelerations . There may be some difference in the 
case of unstable convection in that the critical Rayleigh 
number may be exceeded by the higher impulsive accelerations 
and not by the lower level steady accelerations. 

Vehicle rotations , in addition to producing steady low- 
level acceleration described previously, will introduce 
rotational flow in a contained fluid. These flows are intro- 
duced each time there is a change in angular momentum and 
decay by Ekman damping. Damping times depend on the size and 
shape of the container, viscosity of the fluid, and change in 
angular rate. These times typically range from seconds to 
tens of seconds. 

The types of material processing experiments being con- 
sidered for Shuttle can be grouped into four categories; 

(1) contained solidification experiments, (2) quasi- 
containerless experiments, (3) containerless experiments, and 
(4) fluids experiments. 

Contained solidification experiments are characterized 
by relatively small dimensions (typically, a centimeter or 
less) and high thermal gradients. This will result in 
Grashof numbers on the order of unity for 10 g. Generally, 
the emphasis will be on establishing and maintaining a 
diffusion-controlled interface. Therefore, it is important 
to avoid circulating flows that would disturb this diffusion 
layer. Even though the induced flows are low, if the 
accelerations act continuously in the same direction for 
tens of minutes, displacements on the order of the con- 
tainer dimensions can occur. Therefore, a careful analysis 



of this effect is necessary and some thought should be given 
to the possibility of orienting the sample relative to the 
vehicle so that the flows from sustained lowr'level unidir- 
rectional forces will be minimized. Similar problems could 
also arise from uncompensated transient accelerations. 
Orienting the experiment in a thermally stabilizing configu'- 
ration along the Shuttle + z»-axis (up relative to the payload 
bay) and flying in a gravity gradient mode with the belly in 
the ram position would minimize flow effects from both aeroi-" 
dynamic drag and uncompensated thruster firings. 

Compensated transients arising from internal forces will 
probably not be significant because of the small net dis- 
placements. Rotational changes, however, should be avoided 
during such experiments to prevent circulating flows. 

Quasi-containerless experiments such as floating zones 
or pedestal melts are characterized by a fluid (usually a 
melt) supported at one or more points by a solid connected 
rigidly to the vehicle. Such systems generally have the 
same dimensions of contained solidification experiments and 
are subject to the flows described previously. However, the 
presence of a free surface makes such systems also subject to 
disturbances from compensated transient accelerations which 
will cause oscillations of the fluid relative to the support. 

Container less experiments must have a weak suspension 
system using acoustic, electromagnetic, or electrostatic 
forces to position the sample. Such a suspension system must 
be designed to compensate for residual accelerations such as 
drag, Keplerian effects, vehicle rotations, and various 
transient accelerations. In such experiments the emphasis is 
on maintaining the sample position without physical contact 
with the container walls. Therefore, the low-level steady 
forces are of less consequence, but the transients become 
important, especially high-level transients that can arise 
from internal disturbances, e.g., crew activity. 

Generally, the suspension can be characterized as an 
energy well described by ip (x) , The maximum force that can be 
exerted by the suspension is given by -i|^ 'max(t) and the 
system must be designed such that this is greater in magni- 
tude than the largest expected sustained force. Transient 
forces greater than lil^'maxC^^) I can be tolerated provided the 
total impulse does not impart sufficient kinetic energy to 
the sample to enable it to climb out of the potential well. 
Therefore, it is the impulse, not the peak acceleration, that 
is of concern to the success of containerless processing 
experiments , 
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Fluid experiments encompass a wide variety of experi'- 
ments that are generally characterized by larger dimensions 
but much lower thermal gradients than encountered in con?- 
tained solidification experiments. Such experiments include 
crystal growth from aqueous solution or vapor, electro-' 
deposition, electrophoresis, suspension polymerization, etc. 
Often the emphasis is on maintaining suspensions of multi-' 
phase systems as well as on controlling density-driven con- 
vection. In many cases it is desired to observe the nature 
of convective flows and correlate their behavior with 
vehicular accelerations and experimental effects. 

Such experiments could be seriously affected by low- 
level steady accelerations, uncompensated transient 
accelerations, and vehicular rotation, as in the case of the 
contained solidification experiments. Compensated transients 
are not expected to be very significant. Although sedimentation 
of phases with different densities will be slowed by six 
orders of magnitude from earth-based experiments, it may be 
desirable to avoid long periods of low-level unidirectional 
acceleration such as would be produced by atmospheric drag 
in an earth-fixed orientation. 

In summary, there are multiple factors that must be con- 
sidered in the acceleration environment of a space vehicle 
whose importance depends on the type of experiment being 
performed. This is illustrated in the matrix shown in 
Table I . Some control of these factors may be exercised by 
the location and orientation of the experiment relative to 
Shuttle and by the orbit and vehicle attitude chosen for the 
mission. The effects of the various residual accelerations 
can have serious consequences to the control of the experi- 
ment and must be factored into the design and operation of 
the apparatus. 
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CONVECTION IN FLUIDS AT REDUCED GRAVITY 


by 

Simon Ostrach 

Case Western Reserve University 
Cleveland, Ohio 


At the start of the materials processing in space pro- 
gram it was usually considered that at "zero gravity" there 
would be no fluid motions due to buoyancy. Since such flows, 
called natural convection, were considered to be deleterious, 
there was great enthusiasm for processing materials at "zero 
gravity." Firstly, I want to say that the gravitational 
environment in a spacecraft is of the order of 10“4 g to 
10“^ g, and this is not zero! Therefore, I sincerely hope 
that the term "zero gravity" will be eliminated from our 
thoughts and publications. It should, more properly, be 
replaced by micro-gravity, reduced-gravity, or low-gravity. 
Secondly, it is important to note that natural convection is 
not always harmful and, therefore, to be avoided. In some 
situations it may be desirable to have fluid flows in space 
processes, e.g. , to stir the fluid phase for mixing and 
cooling or to help maintain concentration gradients. In any 
event, it is important to know the extent and nature of con- 
vection in space and the factors on which it depends, in 
order either to minimize the effects of convection, or to 
utilize the convection to advantage. 

Unfortunately, convection, even in a normal environment, 
is one of the more complex fluid phenomena, so it is diffi- 
cult to make specific predictions. The space environment 
further complicates the situation because there are a variety 
of nongravity forces that can also induce fluid motion. Such 
nongravity driving forces which are usually suppressed by 
gravity include surface or interfacial tension, thermal 
volume expansions, g- jitter — and that's another term 
that I'd like to see done away with because I think it leads 
to some misconceptions — and magnetic and electric fields. 
Furthermore, gravity-induced convection can still be appre- 
ciable even at 10 "^ g and lower under certain conditions. In 
order to assess properly the effects of natural convection in 
a reduced gravity environment, it is essential to comprehend 
fully what we know about it at the present time. Thus, I 
will outline for you today some of the key aspects. 

Most of the configurations related to materials proces- 
sing are ones in which the fluid is confined by rigid 
boundaries. Internal convection problems are considerably 
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more complex than external ones because the fluid boundary 
layer and the core are closely coupled. This constitutes the 
main source of difficulty in predicting the resulting flow 
and transport. More than one core configuration, that is, 
more than one flow, is sometimes possible, for a given set of 
conditions, and which will actually occur for a given set of 
conditions, cannot be determined a priori . Furthermore, the 
entire flow is sensitive to configuration geometry and to the 
imposed thermal boundary conditions. 

To add to the complexity it's important to note that 
there are essentially two basic modes of flow generated by a 
body force. The first is what I term "Conventional convection 
(Figure 1). A flow starts immediately. It may be that the 
flow is sufficiently slow so that there is no heat or mass 
transfer, but this is not an instability mechanism. This is 
an immediate flow and yet it has at times been considered to 
be due to an instability. 

The second basic mode is one where the density gradient 
is parallel to the body force, but opposed to it. This is 
what I call "unstable convection." 

In this situation, we have a fluid in the state of 
unstable equilibrium and the flow will start only when the 
density gradient exceeds a certain critical value. In that 
case there is a transient period, and, finally, a steady flow 
is established which takes on the form of B^nard cells or 
transverse or longitudinal vortex rolls, or some combination 
of them. 

Once this is comprehended, you realize what information 
is needed to assess convection. You have to know both the 
magnitude and direction of acceleration, and certainly, in 
some of the early spacecraft flights, only the magnitude was 
reported— the direction is equally as important. You have to 
know the geometric configuration, you must know the boundary 
conditions, and you must know the material properties. 

As if all of this were not sufficiently difficult to deal 
with, both modes of convection, conventional and unstable, 
can occur simultaneously in a given configuration. 

Further complexities arise in material processing 
because the fluid is not homogeneous and, therefore, we have 
concentration gradients that can cause density gradients in 
addition to the temperature gradients. The density gradients 
due to those two causes can enhance or oppose each other in 
coupled ways and relatively little work has been done on such 
problems. In the unstable configuration where you have a 
temperature and concentration gradient parallel and in 



different directions, you can have any number of situations. 
For excimple, a temperature gradient can be stabilizing and a 
concentration gradient can be destabilizing. The density 
gradient, however, can be stable so you have a gravitation- 
ally stable situation. There are some gravitationally stable 
situations like this, that, nevertheless, do lead to very 
complex flows. These are the kinds of things that one sees 
in some of the ocean problems where salt fountains and salt 
fingers occur. Of course, there are some situations where we 
have gravitationally unstable situations. 

Corresponding to the conventional convection, we can 
have concentration and temperature gradients augmenting each 
other or retarding each other. Relatively little work has 
been done on any of this latter type of problem. There has 
been some work, of an uninteresting type, I think, on these, 
where similarity solutions are found and the results are 
quite predictable. 

So, for some of these reasons that I have given here, 
it' s necessary to study each problem more or less by itself 
if you want to obtain detailed information on the transport 
processes. However, you can get considerable insight into 
the qualitative nature of the problem from dimensional 
analysis. The first step in the understanding of these kinds 
of complex phenomena, therefore, is to obtain the relevant 
dimensionless parameters . These are best determined from the 
basic equations and associated boundary conditions. 

The parameters that occur for purely natural convection 
flows are the Grashof number , which equals the product of the 
density difference, gravity, and the cube of a characteristic 
dimension divided by the density and the square of the kine- 
matic viscosity. This is a measure of the relative buoyancy 
to viscous forces. A Rayleigh niomber appears in the energy 
equation, and is the product of the Grashof and Prandtl 
numbers. The latter parameter involves only the thermo- 
physical properties of the fluid. If the density difference 
is due to temperature gradients alone, then the Grashof 
niimber is ggATd^/v^, where 3 is the volumetric expansion 
coefficient. And if the flow is thermally induced, the 
velocities can be estimated from U = /Gr(v/d) for Gr > 1 and 
U = Gr(v/d) for Gr < 1. If the density gradient is due to 
concentration gradients, you have a solutal Grashof number 
and an associated Rayleigh number that looks like 

Gr = 22*^ and Ra - PrGr 
s ^2 s s 
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In this way an estimate of the conventional natural con- 
vection velocities can be obtained. 

For the thermally unstable case, the critical Rayleigh 
number, which is the nondimen sional critical density gra- 
dient, required for the onset of flow must be determined. 

Many years ago some of my students and I showed that these 
critical Rayleigh numbers are very much affected by geometry. 

I have plotted some of the work that was taken from Edwards 
and Catton at UCLA, which shows what happens in cylinders of 
various cross sections (Figure 2) . This figure is for a 
square cross section/ it really does not make any difference 
if it is circular. It can be seen that the configuration is 
essentially an unbounded horizontal layer, H/w = 0/ the 
critical Rayleigh number is 1708, which is a well-known 
niomber. However, as you begin to confine the fluid you can 
raise the critical Rayleigh niimber up to millions. And if 
you look at a more general configuration, like a rectangular 
parallelepiped, in Figure 3 the effect of both aspect ratios 
can be seen. Again, the critical Rayleigh number can be in 
the millions. Thus, it is my general feeling that the 
unstable mode is not going to be a serious one to deal with 
in a reduced-gravity environment, because the Grashof numbers 
or the Rayleigh numbers will be six to seven decades smaller 
than normal; and since most of the problems, in materials 
processing, at least, are confined problems, the critical 
Rayleigh numbers will be relatively high. 

Now, again, as I pointed out before, both modes are 
possible in a given configuration, and to illustrate this, we 
see in Figure 4 a rectangular parallelepiped with one wall 
hot and one wall cold. Now, if this is oriented at 90 degrees, 
which means it is in a horizontal position with the hot wall 
below, you notice that you get no convection until you come 
to a Rayleigh number of 1708; and then, as you rotate it 
around to the vertical position with one wall hot and one 
wall cold, that is 6=0, you have a steady two-dimensional 
unicell. If you rotate it all the way around so that it is 
horizontal with the hot wall above and the cold wall below, 
the configuration is stably stratified, and, again, you have 
no convection. But, intermediate to that, you have the two 
modes interacting and you can have longitudinal rolls, 
meandering flows, transverse traveling waves, and longi- 
tudinal rolls again. So, the situation is fairly complex. 

Now, let us go back and look again at some of the 
detailed kinds of convection and see if we can have buoyancy- 
induced convection with a background g in a spacecraft. The 
acceleration of gravity on earth is usually at a constant 
value of 980 cm/sec^. a steady and reduced value of acceler- 
ation is inherent in most spacecraft because of atmospheric 



drag, centrifugal force due to vehicle rotation, gravity 
gradients, solar wind and solar pressure. Superimposed on 
this background of a \aniform, but weak gravitational field, 
are temporally variant accelerations that are due to engine 
burns, attitude control maneuvers and onboard vibrations from 
machinery or astronaut movement. Some of the steady-state 
accelerations in several missions are shown in Figure 5. 

It can be seen that the background accelerations are on the 
order of 10“5 up to 10"^. To give you a feeling as to whether 
10“^ g or 10”^g is zero, I have computed a Grashof number for 
rather mundane conditions. Suppose we have a temperature 
difference of about 10°K at a temperature level of 20 ®C and 
a characteristic dimension of 10 cm. At 1-g for a gas, the 
Grashof niimbers would be on the Order of 10°, for liquids on 
the order of 10^, and liquid metals on the order of 109. If 
you reduce these by six decades, they still are not zero. I 
have also computed some representative velocities at 10”5 g, 
and in Fig. 5 you see that for liquids and liquid metals you get 
velocities of fractions of a millimeter per second and for 
gases you can get velocities on the order of a half- 
millimeter/second. Of course, if the conditions are more 
severe than that, you can get all other kinds of results. 

As contrasted to the nonuniformities of the gravita- 
tional field, some of which have to be accepted as part of the 
natural environment in the orbital vehicle, there are tran- 
sient or time-variant perturbations to the gravity field at a 
point. The unsteady variations, which are referred to as 
g- jitter, can arise from spacecraft maneuvers and mechanical 
vibration. As I said before, I think g- jitter is an improper 
name, because what we are talking about now are acceleration 
levels of the spacecraft itself, and that is what 
provides the driving force for the fluid. During spacecraft 
maneuvers again, acceleration magnitudes of 10“4 g and higher 
can be encountered. This effect possibly could be overcome 
if experiments were done during the drift mode. On the other 
hand, g- jitter caused by mechanical vibrations cannot neces- 
sarily be controlled at their sources. Mechanical vibrations 
that are transmitted to the experiment have the same effect 
as a time-variant gravity. They can be caused, for example, 
by astronaut movement, rotating and reciprocating machinery, 
extravehicular activities, and arm and leg movements. 

The earliest study of this phenomenon was done at 
Lockheed, Huntsville, by means of computer solutions. 

Numerous instances were considered (see Figure 6) ; a 
rectangle which is heated from the side, a rectangle which is 
heated from below, and a cylinder heated from below. Various 
modes of jitter were imposed on them, such as sinusoidal and 
absolute value of the sinusoid and linear periodic (see 
Figure 7) . The greatest effect of this mode was obtained for 
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the saw-tooth g- jitter (linear periodic) . In Figure 8 the 
isotherms for this g- jitter model are shown to be very simi- 
lar to the isotherms for the model where there is a constant g 
at the same mean level. The streamline patterns are shown in 
Figure 9. This was the first indication of these kinds of 
things, but I do not know how definitive these results are 
because they were obtained merely from large-scale computer 
calculations. A few years ago my group started both theo- 
retical and experimental studies of g-jitter. The dimension- 
less parameters for this problem were found to be as given in 
Figure 10. The basic configuration considered was a rectangu 
lar enclosure which had a hot wall and a cold wall to 
stabilize it convectively so that there would be no motion 
other than that due to the g-jitter. The g-jitter was 
imposed normal to the temperature gradient. The unsteady 
effects are then determined by L^/a. Another parameter is 
the aspect ratio. The next parameter is the volximetric 
expansion coefficient times the temperature difference, 
this is the relative imposed temperature difference, and 
then there is a Reynolds number based on velocity which is 
to X a dimension. 

Now, typical values of the parameters for the given con- 
ditions are also shown in Figure 10. Our analysis indicates 
that if the jitter was of the form fjj = AF(t), where F(t) 
was a random function, but had a zero time mean, like a 
sinusoid, in other words, if the imposed acceleration essen- 
tially had as much on the plus side as on the negative side, 
this would correspond to having a fixed experiment in a space 
craft and having nothing that moves it or moves the space- 
craft at the same time, it was then found that there would be 
no motion ... no significant motion. All that would occur, 
at most, is a sort of random motion with no net translation 
of the particles and no influence on the transport. However, 

if you had a nonzero mean type of acceleration like 

fx = K + AF (t) , where K is a time mean, then, in fact, you get 

an equivalent flow like a natural convection flow where K is 

the essential driving force. There are other aspects of the 
g-jitter problem that were investigated. In the case just 
discussed, the g-jitter acts across a temperature gradient, 
and this corresponds to the conventional convection. If, 
however, the g-jitter was along the temperature gradient, the 
unstable mode of convection would be expected. It appears 
from the work of Greshko and Sani that for that mode, the 
oscillations of the accelerations tend to stabilize the situ- 
ation. There are also other issues involved, such as the 
transient nature of how long does it take for the g-jitter to 
stabilize . 
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Experiments on the conventional mode of g- jitter were 
performed in the Lewis drop tower; the conditions (shown in 
Figure 11) exceeded the values shown in Figure 10* _ All of 
the tests indicated no motions at all. The dimensionless 
parameters were increased in an attempt to get a motion and 
it was at that point that the analysis indicated to us that 
no motion could" be obtained with non-zero mean jitter. 

We just have run a couple of tests where the jitter has 
a non-zero mean, and we did pick up some of the oscillatory 
motion. The results indicate that astronaut movements, like 
sneezing and coughing, which impose jitter with zero 
time means v/ould induce, at most, an oscillatory 
motion which is quite small. There would be no net fluid 
translation and no net heat transfer change. Spacecraft 
maneuvers and astronauts walking along the periphery of the 
spacecraft or violent exercises will lead to non- zero mean 
jitter, and then motion and transfer at the mean level will be 
obtained. An interesting sidelight on this was obtained in 
one of the experiments in which all the bubbles were not 
removed from the liquid. The bubbles caused a very large 
streaming motion, which implies that if foreign particles are 
in the fluid, they could induce rather large scale motions. 

Another possible driving force for fluid motions is 
surface tension. Surface tension on the free surface of a 
liquid may, under some conditions, considerably affect liquid 
motion. The presence of an interface between two fluid faces 
can influence the motion of fluids when either the interface 
has a finite curvature which is different from the equi- 
librium curvature, or when an interfacial tension varies from 
point to point. In both cases, forces appear in the inter- 
facial region that can affect or generate the fluid motions. 

It has generally been believed, and argued that on earth, 
this Bond number (Figure 12) is what indicates the relative 
importance of gravity to surface tension. Computations made 
on this basis indicate that unless the dimension is very, 
very small, the surface tension forces are suppressed. This, 
however, is what I call a static Bond number; in other words, 
it is one which occurs under isothermal conditions, and I 
will show shortly that it is not really a relevant one. The 
relative magnitude of surface tension and buoyancy forces has 
been considered to be given by the ratio of the Marangoni and 
Rayleigh numbers . It will be shortly shown that this is not 
correct either. 

It is important to recognize, and this seems to be very 
much confused, even though a number of years ago Skip Scriven 
pointed this out, that there are two modes of surface tension- 
driven flows just as there are two modes of buoyancy-induced 
flows. If the gradient of concentration or temperature is 



along the interface, then you have a flow immediately 
establishing itself. If the gradient is normal to the inter- 
face, you have an instability, just like a Benard insta- 
bility, now it is the Marangoni instability (see Figure 12) . 

It is very important that this be understood. Now, the 
greatest effort has been given to the unstable type of flow, 
but for materials processing experiments, in fact, the other 
kind of flow is considerably more important. 

I was rather amazed when I was looking at this field 
that the dimensionless parameters had not been formally 
derived and so I set out to determine them. 

Basically, consider a simple configuration of a rec- 
tangular container filled with a fluid with a free surface 
along which a temperature gradient is imposed. Regular 
normalization of variables yields the familiar parameters 
shown in Figure 13. However, this is all in terms of a 
reference velocity, which is, as yet, unknown. Most of the 
people who have looked at these kinds of problems used for 
the reference velocity something like the kinematic viscosity 
over a length or the thermal diffusivity over a length and 
this is not the correct one. The correct one is obtained 
from the free surface boundary conditions . The boundary con- 
dition that must be satisfied (Figure 14) is that the shear 
stress of the wall must be balanced by the surface tension 
gradient. For the configuration under consideration if the 
Reynolds number times the square of the aspect ratio is small, 
we have a viscous -dominated situation, so the characteristic 
dimension, the proper one, is the height of the fluid layer/ 
then the proper velocity is shown at the top of Figure 14. 

This is interesting in that it is given in terms of the 
physically imposed conditions and the thermophysical proper- 
ties of the fluid, viz., the surface tension gradient, the 
imposed temperature difference, the viscosity, and an aspect 
ratio. As a consequence, the important parameter for surface 
tension-driven flows, is what I call a surface-tension 
Reynolds number, which is defined in Figure 14. 

Now for the other case of a boundary-layer situation, 
the boundary layer thickness must be used as the reference 
length and the reference velocity comes out as shown in the 
lower part of Figure 14. Again, the surface tension Reynolds 
niimber is shown thereon. This situation is completely analo- 
gous to what happens in natural convection flows and so if 
you look at the resulting equations shown in Figure 15 for 
the viscous-dominated situation, you notice that the 
Marangoni number appears only in the energy equation. Thus, 
it is not the fundamental parameter that it has been made out 
to be. It is like a Pec let number, which is important for 
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heat transport. For the boundary layer situation (lower 
equations) the Marangoni number does not even appear 
explicitly. 

There are other characteristic lengths for other kinds 
of problems that one has to be careful about, but for this 
kind of geometry, this is what they come out to be. Also, 
very interestingly, you notice that the buoyancy effect is 
given by the parameter shown at the lower left of Figure 15 
and it is now a modified Bond number. For the boundary layer 
case (lower right of Figure 15) it is a modified Bond nimiber 
of a different kind. Once this was recognized, it turns out 
that we have done an analysis and we know that we can now 
look at surface tension flows on earth, because what happens 
is that the characteristic length is determined by the 
imposed buoyancy field and you can properly tailor the flow 
field by heating from above in a stabilizing situation. In 
this way, surface tension flows in fish tanks were observed 
for rather large periods of time. 

One important result obtained from the determination of 
the proper dimensionless parameters is that the surface 
tension Reynolds number in contra-distinction to the con- 
ventional Reynolds number does not depend on the dynamics, 
but depends primarily on the imposed conditions and on the 
fluids. Thus, the problems can be categorized according to 
the types of fluids, and the qualitative nature of the flows 
to be expected in different kinds of fluids can be obtained 
a priori (see Figure 16) . 

Brief mention will now be made of other possible types 
of convection in a reduced-gravity environment (see Figure 
17) . Those that are dependent on the acceleration field are 
listed there. Independent of gravity, thermal volxime expan- 
sions can act as a driving force for convection. After some 
of these earlier Apollo XIV and XVII heat flow and convection 
experiments we thought we had some evidence of this type of 
convection, but I do not think it turned out to be that way. 
In other words, if you rapidly heat a confined fluid, like a 
gas, you have rapid local expansions, which, in turn, generate 
pressure waves. These pressure waves produce a convective 
motion which can greatly increase heat transfer relative to 
conduction and can also cause mass and chemical species 
transport. This provides a mechanism for enhancing or sup- 
pressing convective motions at low gravity by controlling the 
heating rate. 

Phase change convection is another type. Shrinkage 
usually occurs during solidification because the density of 
the solid is usually higher than that of the liquid and, 
hence, occupies a small volume. This volume reduction 
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results in a flow of liquid toward the solidifying interface. 
Such flows are also generated by non-equiraolar reactions that 
occur in vapor deposition crystal growth. Pressure pulses 
can also accompany these kinds of flows and^ therefore, 
influence the nature of the material. Relatively little work 
has been done on this type of convection, and I think it is a 
rather important one that ought to be considered. 

Thermosolutal convection, of course, is one which occurs 
when we have concentration gradients causing the flow, and I 
have already talked about that. It is important to \ander- 
stand that electric and magnetic fields induce body forces, 
just like gravity, so that they can generate convection and 
phase separation phenomena similar to those in the gravi- 
tational field. Whereas, in a gravitational field, convec- 
tion is generated by density differences, in an electrical 
field the flow is generated by differences in electrical con- 
ductivity and by differences in susceptibility in a magnetic 
field. The electric conductivity and magnetic susceptibility 
are temperature dependent, so that temperature differences 
are usually required to obtain these flows. Again, both con- 
ventional and unstable types of convections are possible with 
these kinds of fields, and, of course, in a given problem, you 
can have combined or coupled convections due to any one of 
these . 
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INFORMATION NEEDED TO ASSESS CONVECTION 

(A) BOTH MAGNITUDE AND Dl RECTION OF ACCELERATIONS 

(B) GEOMETRIC CONFIGURATION 

(C) IMPOSED BOUNDARY CONDITIONS 

(D) MATERIAL PROPERTIES 

FIGURE 1 
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STEADY REDUCED-GRAVITY ACCELERATIONS 


SOURCE 

TYPE OF MISSION | 

LOW-ALTITUDE 
EARTH ORBIT 

low-altitude 

LUNAR ORBIT 

INTERPLANETARY 

TRAJECTORY 

ATMOSPHERIC DRAG 

5x10-5g & LOWER 

NONE 

NONE 

CENTRIPETAL FORCE 

1 X 10*®g 

3 X IC'^g 

3x10'^'*g 

GRAVITY GRADIENT 

3 X 10‘® g/cm 

2 X 10“® 9/cm 

4x10'^^ g/cm 

VENTING THRUST 

ALL 10'^g TO 10 ‘®q, 

NEGLECTING DRAG 

VEHICLE THRUST 

ALL 0.1g TO 6g 




STEADY g-LEVELS IN APOLLO MISSIONS 


APOLLO FLIGHT MODE 

TYPICAL 9- LEVEL 

PASSIVE THERMAL CONTROL 

3 X 10*® g 

ATTITUDE HOLD IN TRANSLUNAR OR 
TRANSEARTH ORBIT 

7x10®g 

ATTITUDE HOLD IN LUNAR ORBIT 

5 X TO'^g 


ONE-g GRASHOF NUMBER AT = lO^K, T = 20°C, d = 10cm. 

GAS ^ O (10®) LIQUID ~ O (lO^) LIQUID METAL ~ O (10®) 

VELOCITIES AT 10-5g 

LIQUID AND LIQUID METALS - O (0.13mm/sec.) GAS ~ O (0.5mm/sec.) 

FIGURE 5 
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FIGo 6, CONFIGURATIONS USED IN COMPUTER MODELS 
OF g-JlHER CONVECTION. 
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MODEL 1 
(SINUSOIDAL) 



(ABSOLUTE VALUE SINUSOID) 


Ah—. 


7m/ 


g = MEAN GRAVITY LEVEL 
A = AMPLITUDE OF JITTER 
X = PERIOD OF JITTER. 


t = TlME. 


X 


MODELS 

(LINEAR PERIODIC) 


Fia 7o G-JITTER MODELS USED FOR PARAMETRIC SENSITIVITY STUDY« 
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HOT WALL 



X HOT WAL L X 


a. RECTANGLE (SIDE) b. RECTANGLE (BELOW) 

CONSTANT g 

g- JITTER 

CONDUCTION ONLY 



HOT WALL 

c. CYLINDER (BELOW) 


FIGo a ISOTHERMS AT t - 3 MINUTES FOR THREE CONFIGURATIONS 
SHOWING EFFECT OF q-JIHER MODEL 
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g- JITTER 


PARAMETERS: - 


UNSTEADY 

_ ^ _ _ ^ ASPECT RATIO 

a INERTIA L 2 


ojU2 inertia 

- — L_ =Re ~ - : 

VISCOUS 




g-JITTER 

02 >01 


L-| = 10cm, Li/i _2 = 5, 02 - 0i = 30°C a 0.04 cm. 
OJ/27T « 0.5 RAD/SEC. A = aco^ ~ 10’3 + lO’Sg 

a/L-] = 4 X 10’3, j3 (02 - 0>i) = 10‘2, Re = 4 x 10^ 

Pr = 6 
fx = AF (t) 

•-RANDOM, ZERO TIME MEAN VALUE 
fx = K + AF (t) NONE - ZERO MEAN 
^TIMEMEAN 

FIGURE 10 
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EXPERIMENTS: 

AT = QO - 43°G CO = 10-40 CPS 

a =0.2 cm a/Li = 0.15 

e = /3 AT = 0-1.12 X 10-2 Re = 1.27 TO 5.0 x 106 


RESULTS: 

ASTRONAUT MOVEMENTS ~ ZERO TIME MEAN 
UNICELLULAR OSCILLATORY RATIO O (eaco) 
NO NET FLUID TRANSLATION 
NO HEAT TRANSFER CHANGE 

SPACECRAFT MANEUVERS - NON-ZERO MEAN 
MOTION AND TRANSFER AT MEAN LEVEL 

BUBBLES CAUSED STREAMING 


FIGURE 11 




CONVENTIONAL FLOW 

UNSTABLE FLOW 

SURFACE TENSION 

GRADIENT ALONG INTERFACE 

GRADIENT NORMAL TO INTERFACE 

BUOYANCY 

GRA0I ANT NORMAL TO B. F. 

GRADIENT PARALLEL BUT OPPOSITE TO B. F. 


BOND NUMBER: 


pgd2 GRAVITY 

a SURFACE TENSION 


p - FLUID DENSITY 
g ^ GRAVITATIONAL ACCELERATION 
d - LENGTH 
a - SURFACE TENSION 


FIGURE 12 


90 













FREE SURFACE 


X,U 

X = X/L, y = Y/h, U = U/Ur, v = V/Vr (h/L) 
p = P/pUR2,r = (T-Tc)/(T„-Tc) 
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FIGURE 13 
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FOR Rff a2 « 1 


0 = a2 


92u + 

92u 


9x2 

9y2 

9x 

92y ^ 

92y 

_ GrA ^ 

9x2 

9y2 

Ra 

11 

> 

a2 

92t 92 T 


- (—) 
a 2 \ 3y y 


\ 9x 3yy 
FOR Ra a2 » 1 


9x2 9y2 


u ^ +v ^ ^ + 9^ 

9x 9y \ Raj 3x2 ay2 

u 9v 9v =/A\2/3 32 v ^ a2v 


9y \ Ro / 9x2 ay^ 


9t , 9t 1 

u + V = 

9x 9y Pr 


A 92t ^ 92t 

R„ / 9x2 ay2 


Ma = Marangoni number = Pr Ra 

For Rff a2 « 1 
Gr A _ p|3gL2 _ -g^ 
ba/bl 


For RoA2 » 1 


FIGURE 15 


93 



FOR: (Tw-Tc) = 5QOC, h = L = 10 cm. 




Gr 

Ma 

Bo 

Bo/Ra^^® 

OILS 

o 

1 

6 

102 _ iq7 

102 - 106 

103 

10 

GLASS 

lO-l 

10 

102 

102 

— 

WATER 

106 

io» 

109 

— 

10-2 

LIQUID 

METALS 

105 - 106 

108 - lOlO 

103 - 105 

— 

10-2 


FIGURE 16 


94 


i 

























TYPES OF REDUCED-GRAVITY CONVECTION 

ACCELERATION FIELD DEPENDENT 

THERMAL CONVECTION 
SOLUTAL CONVECTION 
g - JITTER CONVECTION 
THERMOSOLUTAL CONVECTION 


INDEPENDENT OF ACCELERATION FIELD 


SURFACE TENSION CONVECTION (MARANGONi; 
THERMOACOUSTIC CONVECTION 
PHASE CHANGE CONVECTION 

CONVECTION DUE TO ELECTRIC AND MAGNETIC FIELDS 


COMBINED OR COUPLED CONVECTION 

FIGURE 17 


95 


i 



EFFECTS OF SPACECRAFT MOTIONS ON 


FLUIDS EXPERIMENTS 
by 

Roger F . Cans 

Department of Mechanical and Aerospace Engineering 
The University of Rochester 


ABSTRACT 


The equations of motion governing an incompressible 
fluid contained in an orbiting laboratory are examined to 
isolate various "fictitious forces" and their relative influ- 
ence on the fluid. The forces are divided into those arising 
from the orbital motions and those arising from small local 
motions of the spacecraft about its center of mass. The 
latter dominate the nonrotating experiments. Both are 
important for rotating experiments. 

A brief discussion of the onset of time-dependence and 
violent instability in earth-based rotating and processing 
systems is given. 


The dynamics of a fluid in an orbiting spacecraft, 
whether stratified or not, whether contained or not, will be 
affected by all the motions of the spacecraft. I address 
contained unstratified fluids below. 

Experimental data will be measured in a spacecraft-fixed 
noninertial frame. Thus measurements must be interpreted to 
take into account "fictitious forces." 

A spacecraft-fixed coordinate system is shown in Figure 
1. Three orbital modes, gravity gradient, z-local vertical, 
and inertial hold, are shown in Figure 2. 

Both orbital and local (thruster burns, etc.) frame 
accelerations are important. Because they have different 
time scales, I will separate them by making two coordinate 
transformations . 

First, I will transform to the "orbital system," in 
which the center of mass of the spacecraft (CM) is fixed. 

This introduces a fictitious acceleration 
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( 1 ) 


f = -qx (fixR) - R - 2ftxv - a 

where is the instantaneous orbital angular velocity, R is a 
position vector originating at the orbital focus, a dot 
denotes differentiation with respect to time, v is the local 
velocity and a is the apparent acceleration of the focus. 

The quantity QxR is the secular deceleration of the 
orbital motion and is equivalent to the aerodynamic drag 
deceleration. The quantity a is a radial acceleration, and 
is negligible compared to the gravity gradient forces. 

The Coriolis acceleration, 2£xv, exceeds both drag and 
gravity gradient values if lv| is larger than a few mm/sec. 

The total apparent acceleration is then approximately 


-Qx(fixR) - R - 2fixv + a_j^; =-ftxR , (2) 


If I put R = Rq + £/ where Rq is the orbital radius of r a posi- 
tion vector from the CM, then the equations of HeCM and motion 
in terms of the velocity v measured relative to the orbital 
system are 


3 2 

— + 2QXV +V{P/p - I - 5 ^ 



= vV^v 


V*v = 0 • 


(3) 


For a spacecraft in a perfect gravity-gradient or a 
z-local vertical orbit, only the Coriolis force appears. Its 
magnitude can be estimated in terms of U, T and L, which denote 
characteristic velocity, time and length scales. Figure 3 
shows U, T space for L = 30 cm and = 90 min. 

To examine disturbances (attitude corrections, etc.) to 
the perfect orbit, I make a second transformation to a system 
rotating about the CM at o) (t) , which is changing in magnitude 
and direction. The inertial velocity v, the orbital frame 
velocity u and the local velocity £ are related by 

V = u + ~ 3. 
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The pseudoforces associated with this transfomation, 
combined with those from the first transformation, lead to 
the full equations governing £: 

2 

q_l_ + + 2 (g+w) X 2^ + Vtt = V? 2 + ^ 

V*2 = 0 . (5) 

Here F is the perturbation introduced by the experimenter and 
(fixw) X r is called the Poincar4 force. 

For a perfect inertial hold orbit w = ~0.f w = 0, and all 
the effects of rotation not reducible to potential form 
vanish. 

The disturbance term wxr is probably dominated by crew 
with magnitude ~0. 1 cm/sec^ at frequencies near 1 Hz. Cor- 
rections to the motions induced are made by firing vernier 
thrusters at intervals measured in minutes for times less 
than one second. The power spectrum is presumably dominated 
by these frequency bands. 


The Poincare force has a similar power spectrum. The 
disturbance w can be estimated from oofL =0.1 cm/sec^, where f 
is the frequency of the disturbances. Thus the Poincare 
acceleration is ~10~4 cm/sec^ and will generally be negli- 
gible. 

The experimenter then needs to disturb his system in 
ways equivalent to an acceleration of 1 cm/sec^ or greater to 
perform experiments at frequencies coinparable to these. For 
quasisteady experiments smaller disturbances can be used in 
conjunction with data averaging. However, account must be 
taken of "rectified" flows arising from nonlinear interaction 
in the term. 

If the experiment rotates, the equations of motion can 
be deduced by direct substitution. Let £ = w + fxr where y 
is the experimental rotation rate and w the velocity relative 
to the experimental container. 

For any reasonable rotation rate, the terms ux^ and Oxy 
will dominate Qxu. Since the direction of w cannot be 
controlled, it would be prudent to arrange y to be parallel 
to fl. The remaining Poincar^ force, (£xy) x r will now domi- 
nate the wxr term, and the appropriate approximate equations 
are 
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( 6 ) 


Y 



w + w»Vw + 2 y><w +Vir 


r,2 

vV w 


(wxy) X r + F; 

V«w = 0 ' 


Steady precession, which provides a steady Poincar^ 
force, has been examined in some detail [1]. For values of 
the ratio of precession rate Q to rotation rate = Rp) 

and the Ekman number, E = v/wL^, where v,L denote kinematxc 
viscosity and container length scale, small compared to 
unity, the first effect of precession is to tilt the apparent 
rotation axis by an amount Ofe(Rp) . For perfect spheres , and 
cylinders of certain specific aspect ratios, the tilt ampli- 
tude isO'^RpE"^). 

At tilt amplitudes of a few parts in 100, time depend- 
ence appears. The nature of the time dependence depends on 
the container shape. Further increase in Rp leads to violent 
instability in both spheres and cylinders. 


These laboratory examples are at best crude guides to 
what one might expect in space. No studies have been done 
for even a simple time -dependent situation, let alone the 
sort of pitching one might expect in an orbital situation. 
Such studies would be useful as better guides to space 
phenomena. Since the disturbances are apparent in a con- 
tained incompressible fluid, there is opportunity to study 
them on earth, where the conditions of relative acceleration 
can be controlled. 
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TRANSIENT THERMAL CONVECTION IN MICROGRAVITY 


by 

Robert F. Dressier 

Materials Processing in Spaoe Division 
NASA Headquarters 
Washington , D . C . 


ABSTRACT 


The unsteady two-dimensional thermal convection in a 
cylinder due to a transient acceleration is first solved for 
a step-function excitation. From this, the solution is then 
obtained for an arbitrary time -dependent acceleration. The 
solutions are valid for sufficiently low Rayleigh numbers 
and, therefore, relevant to microgravity fields. As an 
example, two graphs are presented for the convection result- 
ing from the movement of an astronaut inside the Shuttle. 

One graph is for a low value of a2/v, the other for an upper 
value, where a is the diameter and v the viscosity. An 
approximate solution is then obtained for the three- 
dimensional unsteady convection in the interior of a cube, 
valid when the Grashof number is sufficiently small. This, 
therefore, gives approximate flow velocities for an arbitrary 
three-dimensional orientation of the acceleration vector 
relative to the hot and cold faces of the cube. The analysis 
can also be applied to obtain any other convective flows such 
as those caused by g- jitter or variable rotation of the 
Shuttle. 


The Materials Processing in Space program is using the 
weightless environment of space for studies that will require 
a fundamental understanding of the response of fluids to 
various types of inertial acceleration excitations. Most of 
these processes have associated temperature and concentration 
gradients (i.e., density gradients) which are necessary to 
control the processes. A better understanding of these fluid 
effects will assist not only materials processing, but also 
help us to specify in advance the g-level constraints which 
the program needs to impose on the space mission. 

The unsteady two-dimensional thermal convection in a 
cylinder with a linear transverse temperature gradient is 
first solved for a step-function acceleration based upon the 
Navier-Stokes equations. From this basic solution, the 
transient flow can then be obtained for an arbitrary time- 
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dependent acceleration. The solutions are valid for a suffi- 
ciently low Rayleigh number, and therefore are relevant to 
microgravity fields important to the behavior of fluids in 
space. 

Figure 1 illustrates the replacement of the convection 
model for a square cross-section by a cylinder, with little 
error for the maximum velocity, because the maximum velocity 
in both cases occurs far from the outer wall. Our basic 
step-function solution approaches a steady-state flow for 
large time which is identical with the steady solution 
obtained by Weinbaum [1] for a cylinder, and almost identical 
with the steady solution obtained by Batchelor [2] for a 
square cross-section. Our step-function transient is 
expressed in terms of the Bessel functions Jq, J]^, J 2 , and 

Js- 


Figure 2 shows for our step-function transient solution 
the velocity distribution (in dimensionless form) versus the 
cylinder radius for various values of time. The limiting 
profile for infinite time is the steady solution identical 
with the Weinbaxim [1] steady solution. The figure shows the 
effect of the boundary layer moving inward through the flow 
with increasing time. 

In all the following examples, the direction of the 
acceleration vector has been chosen to produce maximum two- 
dimensional convection. Figure 3 illustrates an application 
of our general analysis to the specific example of a 200 lb 
astronaut moving across the Shuttle by pushing against one 
wall to start, and braking against the opposite wall to stop. 
The v®(r,t) is our basic step-function solution, and this 
inertial problem is solved by the appropriate linear combi- 
nation of the basic solution with various phase lags. 

All velocities in our solutions are directly propor- 
tional to the magnitude of the imposed acceleration and to 
the percentage density variation across the cylinder. They 
are more complicated functions of the basic parameter a^/v, 
where a is the cylinder radius and v is the viscosity. 

Figure 4 shows the maximxim velocity, which is approxi- 
mately located at r/a = 0.6 resulting from the astronaut 
motion, for two extreme values of a^/v at 2500 and at 100, 
for water. The former case has a = 5 cm with a half-life for 
decay of 170 seconds. The latter case for a = 1 cm has a 
half-life of 7 seconds. In order to show specific values, we 
have taken a 4% density variation here. 0 is the angle 
through which a particle rotates in its circular path. In 
both cases of a^/v, peak velocities are the same at 
12 microns/sec, but 0 approaches a final value of 0.06® for 
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the large cylinder. It peaks at 0.28® for the small cylinder, 
then returns to a smaller value for large t. These quanti- 
tative results suggest that the astronaut motion in some 
cases may not significantly interfere with fluid-crystal 
experiments in Spacelab. 

By contrast. Figure 5 shows peak velocity and angular 
displacement of the fluid when the Shuttle makes a 90® roll 
by firing its main thrusters initially for 0.7 seconds, 
rolling about 90 seconds at l®/sec, then firing reverse 
thrusters to stop the roll. The same extreme values of a^/v 
are used, with the same 4% density variation, and ass^aming 
the fluid cell to be located one meter from the axis of 
rotation of the Shuttle. The upper two graphs give peak 
velocities and the corresponding displacements due to the 
centrifugal force. For the large a^/v, the velocity peaks at 
0.17 mm/sec, and in both cases 0 reaches a large final value. 
The velocity for the small a^/v has attained almost the 
steady-state value after the first 20 seconds. 

The bottom two graphs show the effect of angular 
acceleration. Here the excitation consists of two short 
pulses, but having a much larger magnitude than the centrifu- 
gal acceleration. The peak velocities are high at 0.12 mm/ 
sec, and the peak 0 values are also high, although 0 
decreases after 90 seconds to small final values. 

These velocities and displacements are all sufficiently 
large to suggest that a roll-maneuver should not be made 
during fluid-crystal experiments in Spacelab. 

When the direction of the acceleration vector is arbi- 
trary, a three-dimensional motion will in general be pro- 
duced. Figtore 6 shows the acceleration vector resolved into 
three components. The (0) component will produce no con- 
vection lanless the critical Rayleigh number is exceeded, and 
only for strict alignment of density gradients and acceler- 
ation, which are unlikely for microgravity excitations. 
Gomponents (1) and (2) would each separately produce a two- 
dimensional convection as shown. If the velocity vector 
fields for these two motions are added together , the sum will 
not satisfy the nonlinear Navier-Stokes equations because of 
the cross-product terms in the acceleration. However, if 
these quadratic terms are small relative to the other terms, 
the superposition can be considered as an approximate solu- 
tion. Figure 6 illustrates the symmetric case where the (1) 
and (2) acceleration components are equal. Calculations of 
errors around a closed trajectory show that the approximate 
three-dimensional flow should be accurate when the Grashof 
number is sufficiently small, as indicated. 
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Both Batchelor [2] and Weinbaum [1] give upper bounds of 
about 2000 on the Rayleigh number for which their steady 
solutions should be accurate. In most of our transient solu- 
tions, the velocities are very much lower (usually only a few 
percent, e.g., for the astronaut-motion example, only 6% for 
a^/v “ 100 and 0.25% for a2/v = 2500) than the corresponding 
steady flow; therefore, the range of Rayleigh number for 
validity can usually be greatly increased. An analogous 
increase should also be possible for the Grashof nxamber since 
this was calculated for steady flows. 
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SPACECRAFT DYNAMIC ENVIRONMENTS 


by 

George H. Fichtl 

NASA/George C. Marshall Spafe Flight Center 


It is useful to present a brief summary of the dynamic 
environment that is anticipated to occur during the Shuttle 
Orbiter missions. The dynamic environment of planned Shuttle 
missions has been examined in a relatively detailed manner 
for the Spacelab Missions 1, 2 , and 3. Accordingly, the 
comments that follow relate to these missions, especially the 
Spacelab 3 Mission since it is the first operational flight 
of Spacelab and is a materials processing, low-gravity 
emphasis mission. However, the ideas presented can be carried 
over to other missions. Zero-gravity conditions in Earth 
orbit cannot be obtained in the Shuttle Orbiter, at least on 
the first three Spacelab missions. However, through careful 
planning, the dynamic environment and its effects on experi- 
ments can be minimized. Furthermore, although the dynamic 
environment of the Shuttle Orbiter is to a large degree 
stochastic, it is possible to predict characteristics of this 
environment so that scientists and technologists can plan 
their experiments and mission managers can plan missions with 
a view toward minimizing the effects of spacecraft dynamics 
on experiments. Characteristics of the dynamic environment 
that might be predicted include typical and "worst case" 
values of vehicle acceleration for the anticipated accelera- 
tion sources, typical number of acceleration events, duration 
times of discrete acceleration events, bandwidth of 
acceleration time history, etc. 

The dynamic environment that will occur during a Shuttle 
Orbiter mission can range over many orders of magnitude (for 
example, 10“^ g to 10~3 g on the Spacelab 3 Mission) and will 
be characterized by a frequency bandwidth ranging from 1/T 
(T is the orbit period) to approximately 10 Hz for rigid body 
accelerations, and from approximately 0.01 Hz to audio fre- 
quencies for the accelerations associated with flexible body 
dynamics, and vibration and acoustic sources. Extensive 
infoinnation is available on the former relative to providing 
scientists and technologists with information about space- 
craft dynamics. This is especially true for Spacelab 
Missions 1, 2, and 3. A unified summary of the latter in the 
context of experiment planning and performance is not avail- 
able. The discussion that follows is directed toward space- 
craft dynamics as related to rigid body dynamics. 
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Table 1 provides a brief summary of the major sources of 
rigid body accelerations and associated values that are 
expected to occur during the Spacelab 3 Mission. The aero- 
dynamic drag and gravity gradient forces create a quasi- 
steady acceleration with characteristic time scale on the 
order of the orbit period. Superimposed on this quasi-steady 
acceleration environment are short-period disturbances with 
time scales on the order of a few seconds or less. The most 
important of these disturbances are the accelerations asso- 
ciated with the vernier control system thrusters and crew 
activity. Water diamps, flash evaporators, and venting of 
experiment gases and liquids will create accelerations in the 
10 “ 7 g to 10”^ g range and thus are small (although not 
necessarily unimportant) compared to the accelerations 
created by the vernier control system and crew activity. The 
accelerations associated with crew motions were predicted 
from force measurement data acquired from the Skylab 
missions [1] . Crew activity has acceleration levels with 
bandwidth ranging from approximately 0.3 Hz to about 2 Hz. 

Vernier control system thrusters are used to maintain 
Orbiter attitude. The duration of a vernier control system 
thruster firing event ranges from 0.08 sec upwards to 0.5 sec 
and depends on the attitude deadband setting. Typically the 
larger the selected attitude deadband, the larger the vernier 
control system thruster firing event time. Thus, at the 
minimum attitude deadband of 0.10 the typical vernier control 
system thruster firing time is 0.08 sec. As the attitude 
deadband is increased, the thruster event firing time will 
increase. The thruster firing rate depends on the vehicle 
attitude and the attitude deadband setting. For a given 
attitude deadband setting, the minimum number of thruster 
firings will occur for the gravity gradient flight mode. To 
minimize the thruster firing rate on the Spacelab 3 Mission, 
a gravity gradient flight mode will be used with the Orbiter 
tail toward the Earth with the wings essentially in the orbit 
plane. On the Spacelab 3 Mission, attitude deadband setting 
will be determined by pallet experiment pointing requirements 
and fluid flow effects that result from spacecraft acceler- 
ations associated with thruster firings. Spacelab 3 Mission 
flight simulations, wherein the effects of random crew motion 
are included, predict 25 to 80 thruster firings per orbit for 
a 1® deadband and 90 to 350 firings per orbit for a 0.1° 
deadband, depending on the crew activity level, so that dead- 
band setting can have a pronounced effect on vernier thruster 
firing rate. The associated number of vernier thruster 
firings per orbit to maintain an inertially fixed attitude 
will be significantly larger than those noted above. 



The number of acceleration disturbances per orbit asso- 
ciated with crew motion is at least a factor of ten greater 
than that associated with the vernier thruster firings. How- 
ever, the spacecraft accelerations associated with crew 
motions do not on average (in time) impart net momentum to 
the fluids, while thruster firings do impart momentiim on 
average. The former is a result of conservation of momentum. 
The reason for the latter situation results from the thruster 
geometry depicted in Figure 1. When the vernier thrusters 
are fired to produce roll, pitch, and yaw motions, both a net 
torque and a net force occur. The net torques can average 
to zero over a sufficiently long period of time, while the 
net forces associated with thruster firings to induce roll 
and pitch modes of motion will not average to zero no matter 
how long the averaging period. This is a direct result of 
the fact that vernier thruster firings associated with pitch 
and roll nations always result in a rectilinear acceleration 
vector directed out of the Shuttle Orbiter bay. Net thruster 
forces associated with yaw motions will average to zero when 
averaged over a sufficiently long period of time. 

As noted above, additional accelerations can occur as a 
result of vibrations and acoustic fields resulting from the 
operation of machinery, etc. A prediction of these environ- 
ments has of yet not been made. However, current constraints 
on acoustic and vibration levels require that the Spacelab 3 
Mission experiment developers design equipment so that 
acoustic and vibration levels of experimental apparatus do 
not exceed the specified levels in the Spacelab Payload 
Accommodation Handbook [2] . 

The Spacelab 3 Mission timeline has been planned so as 
to exclude major vehicle maneuvers. The technical motivation 
for this constraint on maneuvers results from the pronounced 
fluid flows that can occur in a contained fluid when the con- 
tainer walls suddenly undergo a change in rotation. If the 
Orbiter is set into rotation impulsively with angular radian 
rotation rate Afi, then a fluid in a container of characteristic 
size r will experience a flow adjacent to the container walls 
with typical initial velocity rAJ2. This wall boundary flow 
will be transmitted into the interior fluid via viscous 
diffusion of vorticity and Ekman boundary layer effects. For 
a major maneuver with Afl = 2TT/120 rad sec”l and a fluid con- 
tainer with r = 10 cm the initial flow adjacent to the 
container wall, for an impulsive initiation of the maneuver, 
is approximately 0.3 cm sec“l. 

The requirement for constraining the Orbiter attitude to 
one in which the wings are essentially in the orbit plane 
arises from the requirement to preclude the occurrence of 
precessional modes of fluid motion in the spherical 
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convection experiments in the Geophysical Fluid Flow Cell 
(GFFC) [3] as a result of precession of the GFFC angular 
velocity vector about the orbit angular velocity vector. 
Constraining the attitude of the Shuttle Orbiter so the wings 
remain in the orbit plane results in the GFFC angular 
velocity vector being parallel to the orbit angular velocity 
vector. This precludes the occurrence of precessional fluid 
motions in the GFFC. 
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TABLE 1 


ACCELERATION G- LEVELS AND SOURCES* 


Source Disturbance Level 

— q “ _“ o 

Aerodynamic Drag 10 g - 10 g 

Gravity Gradient Force 10" ^g - 10”^ g 


Comments 

Quasi-steady 


Left, Right Flash Evaporators 2 x 10 
Both Flash Evaporators 1 x 10 
Water Dump 4 x 10 
Vernier Thrusters 5 x 10 


- 6 . 

-7 

-6 

-5 


g 


g 

■g 


5 X 


! Dirac Delta 
Function- like 
Disturbances 


Crew Motion (nominal) 10”^g 

_ Z 

Crew Motion (vigorous) 10 ■ g 


Frequency 
■ Bandwidth 
,0.3 to 2 Hz 


*The values in this table are appropriate for 
Spacelab 3 Mission. 
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ORBITAL ENVIRONMENT DYNAMICS 


by 

Robert E. Smith 

NASA/ George C. Marshall Space Flight Center 


Very little was known about the upper atmosphere prior 
to the advent of satellites. During the tracking of the 
earlier satellites, it was found that the upper atmosphere 
was not static but had several variations. Figure 1 shows 
that the atmosphere has a diurnal variation which also varies 
with height, from a factor of three at the lower altitudes to 
a factor of seven at the higher. Figure 2, which is a 
picture of the atmospheric density, shows that the maximum 
density occurs at about 1400 LST, while the minimum density 
occurs slightly before 0400 LST. This maximum density bulge 
follows the sun, so that it is about 2 hours behind the sun 
and at the same latitude; therefore, if the maximum density 
is in the northern hemisphere, then the minimum nighttime 
density would be in the southern hemisphere. In the real 
atmosphere the ratio of the daytime maximum to the nighttime 
minimum density is somewhat a fiinction of solar activity, but 
the atmospheric models, for simplicity, used a constant 20% 
difference between the two values. In the models the con- 
stituents are assumed to be in diffusive equilibrixam above 
about 100 km and all mass densities have been tied to temper- 
ature, so that when the temperature at any altitude is known, 
then there is one density associated with that temperature, 
no matter where it occurs on the globe. The temperature and 
density also are dependent upon the level of solar activity 
and sunspot cycles. Figure 3 shows that the temperature of 
the exosphere varies from around 600®K to on the order of 
2,000®K. In the models all temperatures and densities go to 
a single value around 90 to 100 km. At orbital altitudes the 
models are fairly representative of the real atmosphere; 
however, below 180 km, care should be taken in their use 
because the models aren't very representative of the real 
atmosphere. Variations occur there that are not modeled very 
well because this is a region where there are not too many 
measurements. Single values of temperature, pressure, and 
density at 90 to 100 km are used and the models built up from 
those . 

The upper atmosphere density also has radical changes 
associated with geomagnetic activity. Figure 4 shows that 
the higher the altitude, the more the effect of the geo- 
magnetic disturbance. The ap index is a measure of the 
magnetic activity of a solar storm. In the models the 
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response of the upper atmosphere to a magnetic disturbance, a 
density increase, is Usually global; however, the atmosphere 
itself really responds with a large density increase in the 
auroral zone where the energy from the magnetic disturbance is 
dumped and then a series of waves propagates toward the 
equator. The time delay between a magnetic storm and the 
atmosphere response is very short at the auroral zones; how- 
ever, it takes about 6-1/2 hours for the density waves to 
propagate to the equator. This effect is very short-lived, 
and a spacecraft sees it as a spike in density. 

Figure 5 shows that at 600 km there is about two orders 
of magnitude difference in the density between high solar 
activity and sunspot minimum. 

Some seasonal effects have been found, but they are very 
minor compared to the rest. Changes in atmospheric density 
were correlated with the 2,800 MHz, 10.7 cm, solar flux 
measured at Ottowa, Canada, and also with the geomagnetic 
activity index, although neither of these really affects the 
upper atmosphere. It is the ultraviolet solar radiation that 
does cause the observed changes. Figure 6 shows that these 
indicators are fairly representative of the EUV input and the 
temperature changes which have been observed. 

Recent studies have shown different compositions than 
those portrayed by the older models, while the total mass 
densities have stayed the same. If the constituents of the 
upper atmosphere are important in your application the newer, 
models should be used. 
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DIURNAL DENSITY VARIATION 

THE MAGNITUDE OF THE VARIATION IS INDICATED BY THE VALUES OF TYPICAL MAX TO 
MIN RATIOS NOTED UNDER EACH PEAK. THE PHASE IS INDEPENDENT OF ALTITUDE. 

(ARER JA CCHIA.) 
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TEMPERATURE D I STR I BUTION ABOVE THE THERMOPAUSE WHEN Tq- 1000 K. 
HOURS OF LOCAL TIME, COUNTED FROM MIDNIGHT, ARE MARKED ON THE EQUATOR. 

(AITOFF'S EQUAL-AREA PROJECTION,) 
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EXTREMES OF TEMPERATURE VARIATION OVER THE 11-YEAR SOLAR CYCLE 



4 


Fig. 




ENHANCEMENT OF UPPER ATMOSPHERIC DENSITY DURING A GEOMAGNETIC STORM. 
DATA DEDUCED FROM SIMULTANEOUS OBSERVATIONS OF DRAG ON THE INDICATED 
SATELLITES. THEMAGNETIC INDEX ap IS A MEASURE OFTHEMAGNETIC DISTURBANCE, 
ZERO BEING QUIET TIMES AND 400 BEING A VERY INTENSE STORM. (AFTER JACCHIA.) 
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DAY AND NIGHT DENSITY PROFILES IN 
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COMPARISON OF THE EUV FLUX WITH THE EXOSPHERIC TEMPERATURE, 
THE 2800 Mc/sec SOLAR FLUX AND THE GEOMAGNETIC INDEX 
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Fig. 6 
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CONTRIBUTIONS OF ORBITAL MECHANICS TO 


CONDUCTING EXPERIMENTS IN SPACE 
by 

Larry D. Mullins 

NASA/ George C. Marshall Space Flight Center 


ABSTRACT 


Orbital Mechanics as a discipline is principally con- 
cerned with solving the set of equations (1) for analyzing 
the motion of a satellite under various conditions. This 
activity on the surface may not seem crucial to conducting 
experiments in space, but it provides insights into the way 
in which forces may influence these experiments. More 
directly, for experiments concerned with external "targets,” 
it provides predictions of the satellite's position and 
velocity versus time, enabling extensive preflight planning and 
resulting in optimum use of on-orbit time. 


Orbital Mechanics applied to artificial earth satellite 
motions is concerned primarily with analyzing the motion of a 
point mass under the influence of forces acting on it. Even 
though an earth satellite such as the Space Shuttle may be an 
extended body, it is treated in orbital mechanics as though it 
were a point mass. The satellite's position and velocity or 
state vector is commonly described by a set of six orbital 
elements. For a pure two-body motion (inverse square radial 
force) all of the elements remain constant in time except the 
mean anomaly, which increases at a uniform rate. Additional 
(perturbing) forces cause the orbital elements to change in 
time. This change in time is described by the Lagrange 
Planetary Equations. These are written in equations (1) in 
terms of perturbing forces S,T,N along radial, tangential and 
normal directions, respectively, as illustrated in Figure 1. 
The main problem of orbital raechanics is to describe accu- 
rately the forces acting on the satellite (Figure 1) and then 
to integrate the set of equations (1) . The integrals of 
equations (1) describe the orbital elements as a function of 
time. This task, in itself, may appear on the surface to be 
only indirectly related to the conduction of experiments in 
space . 
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Orbital Mechanics, however, can offer insights into 
planning and analyzing experiments in space by providing 
estimates of the forces acting on the satellite (and experi- 
ment) and describing how those forces affect the orbit of the 
satellite. Many experiments are conducted in space primarily 
because of the (nearly) force-free environment. To be pre- 
cise, the nearly force-free environment is due not to the 
absence of forces on the satellite but to the balancing of 
oppositely directed forces; principally, the balancing of the 
inwardly directed gravitational force and the outwardly 
directed centrifugal force. To further clarify the problem 
one must know if the experiment is rigidly attached to the 
spacecraft, in which case it experiences the same forces (and 
accelerations) as the center of mass of the spacecraft. If 
the experiment is suspended freely from the spacecraft (as, 
for example, a cloud experiment), the most important aspect is 
the relative acceleration between the spacecraft center of 
mass and the suspended experiment. 
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There are four principal types of perturbing forces 
acting on an earth satellite. They are (1) gravitational 
forces, (2) atmospheric drag forces, (3) solar radiation 
pressure, and (4) electromagnetic (Lorentz) forces. These 
forces affect the spacecraft and its contents differently. 

The Gravitational Force 

The gravitational force acts alike on both the space- 
craft and the experiment inside, thus producing no relative 
acceleration between them except possibly for small gravity 
gradient effects. The sources of perturbing gravitational 
forces are the distant celestial bodies (sun, moon and 
planets) and the nonspherical mass distribution in the earth. 
For close earth satellites the latter is by far the largest. 
The forces produced by the celestial bodies can be ignored in 
all but the most sensitive of experiments. The gravitational 
force produced by the earth can be written as the gradient of 
a potential function 


F = - VV 



1 9v , 2 1 iv 

r 86 * r cos 6 d<j) 
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This can be written in terms of S, T, and N by the relations 
r = S ; 6 = cos 6T + sin 6 N ; $ = sin 6 T - cos 0 N, (3) 
as illustrated in Figure 2, where 
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The coefficients (zonal harmonics) and and 

(Tesseral harmonics) can be found to high order in various^ 
publications, e.g. , the Smithsonian Standard Earth (III) 
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(SAO Special Report 353) . The values of the first few coef- 
fieients are 


J 2 = 1.0826x10"^; = -2.546xl0”^j = -1.649xl0~® 

€ 22 = 1.57x10"®; S 2 2 = -.897x10”®; ( 6 ) 

14 1 2 

\x = 3.987012xl0'^^ Ta/s 


Equations (2) through ( 6 ) can be used to ealculate the magni- 
tude of perturbing gravitational forces (those not balanced 
by centrifugal forces) on the satellite. Gravity gradient 
effects are those produced due to the magnitude of the gravi- 
tational force being slightly different at different 
locations in the satellite. 

The principal effects of the gravity perturbations on 
the orbit plane are caused by the first zonal harmonic J 2 . 

It causes a secular regression of the line-of-nodes amounting 
to about 5°/day and short period oseillations in the semi- 
major axis of amplitude up to 7 or 8 kilometers. The J 3 
zonal harmonic causes a long period oscillation in the 
eccentricity of near circular orbits. 

The Atmospheric Drag Force 

The atmospheric drag force is a noneonservative force 
acting generally in a direction opposite the velocity vector. 
This force takes energy out of the satellite orbit, causing 
the orbit gradually to decay. The magnitude is given by 


D 





(m/s^) 


(7) 


where 


3 

p = atmospheric density (kg/m ) 

V = satellite velocity (m/s) 

= coefficient of drag 

(dimensionless number with a value near 2 ) 

2 

A = cross-sectxonal area of satellite (m ) 

M = mass of satellite (kg) 
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The direction of this force, for near circular orbits, can be 
taken as being in the -T direction. The magnitude of the 
force decreases strongly (exponentially) with increasing alti- 
tude because the atmospheric density varies this way. 

Unlike gravity, atmospheric drag exerts a force only on 
those surfaces exposed to the atmosphere. Thus, it would 
exert a force on a satellite but not on a freely suspended 
experiment inside, producing a relative acceleration between 
them of the magnitude given in equation (7) . This would cause 
the freely suspended experiment or sample to drift relative 
to the spacecraft until it contacted some restraining surface 
such as a wall. The magnitude of the drag accelerations on a 
shuttle can vary from as much as 10“4 xa/s^ for a very low 
altitude orbit, say about 250 km, to as little as 10“7 m/s2 
for a very high altitude orbit, say about 600 km. 

Solar Radiation Pressure 


This force is due to sunlight falling on the satellite 
and, like the atmospheric drag force, it acts only on those 
surfaces to which it is exposed, thus producing relative 
accelerations between spacecraft and freely suspended objects 
inside. Its direction is taken to be opposite that of the 
solar vector and it exhibits the unique behavior of "switching 
off and on" when the satellite enters and exits the earth's 
shadow. Its magnitude is given by 





(1+3) E 

£l 

G 



( 8 ) 


where 

3 = the reflectivity (0 < 3 < 1 ) 

O 

c = speed of light ('' 3 x 10 m/s) 

2 

A = area of satellite presented to the sun (m ) 

m = mass of satellite (kg) 

E = "solar constant" ('■1.4x10^ 

^ m^-sec 


For the Space Shuttle (A/m) is on the order of 10“3 m2 /kg, 
making this acceleration on the order of 10”8 m/s2. The 
direction (illustrated in Figure^SJi is taken to be opposite 
the solar vector which in the S,T,N coordinate system is 
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-cos 6* cos a* 
-cos 6* sin a* 
-sin 6* 


(9) 


— 

S,T,N 
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where (fi)z is a rotation about the z-axis through the angle 
Q, etc. 


Electromagnetic Force 

The Lorentz force is that acting on a charged particle 
moving in a magnetic field. Satellites can and do build up 
charges on their surfaces, and the earth about which they are 
orbiting possesses a magnetic field; thus they experience the 
Lorentz force . The magnitude and direction of the resulting 
acceleration is 


a = ^ (V X B) (10) 

m 


where 

q = net charge on satellite 
m = mass of satellite 
V = velocity vector of satellite 
B = geomagnetic field 

The major difficulty in calculating this force is knowing 
what charge is built up on the satellite. In any case this 
acceleration is very small . One source has estimated that it 
is no larger than 10”9 m/s2 in a maximum case. 
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FIGURE 1 orthogonal DIRECTIONS AT THE SATELLITE POSITION: 
S-RADIAL, T-TANGENTIAL, N-NORMAL 


N P. 



FIGURE 2 RELATIONSHIP BETWEEN AND S, T, N 



N. 



FIGURE 3 RELATIONSHIP OF SOLAR VECTOR TO S, T, N 
COORDINATE SYSTEM 
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DYNAMIC BEHAVIOR OF PARTICLES 


IN SPACECRAFT 
by 

Barton S . Perrine 

NASA/George C. Marshall Space Flight Center 


INTRODUCTION 


This presentation is concerned with the behavior ot 
particles relative to a spacecraft frame of reference. The 
discussion was based on work performed by the author in the 
late 1960 's as described in NASA TM X-53643 [1]. 

With the advent of sophisticated scientific space 
missions, the need to understand the behavior of the motion 
of particles relative to a spacecraft frame of reference has 
become of paramount importance. Intuitively, it seems that 
there would be no problem caused by the extremely small 
differential forces acting on the spacecraft arising from 
crew motion, vernier thrusters, small variations in venting 
thrusts, solar pressure and drag decelerations. However, as 
Figure 1 shows, these forces, typical of those acting on a 
spacecraft of differing configuration and orientation, can 
cause large separation distances in a relatively short period 
of time. For this particular figure, the difference in the 
drag deceleration between the two bodies was 0.312 x 10~6 
m/sec2. The technical motivation for the need to understand 
free particle behavior arises from the fact that one may 
require free particles to remain in a specified volume of 
space within the spacecraft. Movement of the particles out- 
side of this volume could result in adverse experimental con- 
ditions, like for example the particle leaving a field of 
view, or striking an experiment boundary. Specific possible 
applications of this knowledge can be cited relative to the 
Drop Dynamics Module scheduled to fly on Spacelab Mission 3 
in 1984, Materials Processing in Space containerless process- 
ing experiments, certain kinds of relativity experiments, as 
well as others. We shall not attempt to present a completed 
or exhaustive theory of particle motion in space, but rather 
we shall outline the kinds of theoretical considerations that 
we must make relative to the definition of experiments. 
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DERIVATION OF THE EQUATIONS OF MOTION 


Let the spacecraft be in a circular orbit of radius, R, 
and angular velocity , oo (Figure 2) . Set up a two-dimensional 
coordinate system in this spacecraft with the positive y-axis 
pointing outward along the radius vector and the x-axis in 
the opposite direction of the inertial velocity vector. Let 
the_position vector in this coordinate system to a second body 
be r. If R is the radius vector from the center of the earth 
to the origin of the relative rotating coordinate system, 
then the differential equation of motion for the center of 
mass of the second body can be written: 

r = a - 2wx r + g - wx [wx (r + R) ] 

where a is the sum of all the differential accelerations 
between the two bodies due to aerodynamics, vernier thruster 
firings/ venting of gases, crew motion, etc., -2o)xr is the 
Coriolis acceleration, g is the gravitational acceleration, 
and the last term is the centrifugal acceleration. The last 
two terms can be combined in the following manner: 

- _ p(r + R) 

"-'IFTiF 

according to the inverse square gravitational law, where u is 
the gravitational constant for the earth. Since o)2 = y/R3^ 

- _ -w^(r + R)R^ 

Ir + R|^ 

Since w is perpendicular to (r + R) , 
wx [wx (r + R) ] = (r + R) 


and 


g - wx [w + (r + R) 3 = (r + R) | 1 - 


Ir + Rp 
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The term | r + R I can be approximated by R + y with an error 
of only about x^/R, where x and y are the components of r in 
Figure 2. Therefore, 


g - cax [wx (r+ R) ] = o)^ (r + R) 



(R + y) 


03^ (r + R) 


1 - 


(1 t |) 


Expanding the second term in the brackets by the binomial 
expansion yields 


g - wxIwxCr + R) ] = (r + R) 


1 - 






where r + R/R^may be approximated by the unit vector in the 
y-direction, j. 


-1 JL 2 ^ 

r = a - 2wx r + 3w yj 

Separating this acceleration into x and y components yields 

X = a + 2wy 

#2 

y = a^ - 2wx + 3w y 

For this report, it will be assumed that the difference 
in drag deceleration acts on the spacecraft in the positive 
x-direction and is a constant, D. The Rand Corporation 
investigated the equations for a drag which is dependent on 
altitude _ [2] . However, this assumption is justified for very 
nearly circular orbits in a spherically symmetric atmosphere 
(with no diurnal bulge) . The governing differential equa- 
tions therefore become 
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X = D + 2wy (1) 

y = -2o)x + 3w^y (2) 

These equations may be solved analytically as shown in 
Reference 1 to give x, y, x, and y as functions of time. 

3 2 

X = x^ + c^t - j Dt - 2b^ sin wt + 2aQ( 1-cos wt) (3) 


y = o 77 ~ “ — t - b^ cos wt + a^ sin wt 
3(0 (0 o o 


X = c^ - 3Dt - 2cob cos (ot + 2coa sin (ot 
o o o 


y = - — - + oob_ sin cot + (oa cos (ot 
(0 o o 


where 


1 • ^ 2d 

(0 -^0 (0 


b = 3y - 2 
o ■'o 


6(oy^ - 3x^ 
o o 


and Xp, yQ, Xq, and y. are the conditions of position and 
velocity at time = 0. 

These equations may be used to describe the motion of 
any body in a near-eircular orbit which remains fairly close 
(i.e., much less than the radius of the orbit) to the origin 
of the rotating relative coordinate system. Also, these 
equations are not restricted to describing the motion of only 
one body. By having a separate set of equations as (3) 
through (6) for each body, the motion of multiple bodies in 
low eccentricity orbits can be described. By subtracting 
their coordinates, the motion of all of the bodies can be 
described relative to one of them. This is a possible method 
of removing the restriction that the spacecraft be exactly in 
a circular orbit. 
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BEHAVIOR OF FREE BODY 


In the following discussion, the body located at the 
origin of the relative coordinate system will be called the 
reference body and the other body the experiment body. 

The equations of motion developed in the previous 
section predict the experiment body to follow a trajectory 
relative to the reference body similar to that shown in 
Figure 3. In this particular example, the motion of the 
experiment body is initiated at the position x = 30 m, 
y = 4 m, with a positive x and y such that its orbital energy 
is greater than that of the reference body. Superimposed 
onto an oscillation in the vertical and horizontal direction, 
the experiment body first experiences a drift to the right 
and then a drift to the left. The oscillation in the 
vertical direction results as it moves between apogee and 
perigee. The oscillation in the horizontal direction is due 
to the fact that the perigee and apogee velocities are 
different. The drift to the right is caused by the semi- 
major axis of the conic of the experiment body being 
initially greater than that of the reference body and thus 
having a longer period. 


CONCLUSIONS 


Significant spatial excursions of particles in space can 
occur relative to the spacecraft frame of reference as a 
result of drag deceleration of the vehicle. These vehicle 
excursions tend to be large as time increases. Thus, if the 
particle is required to remain in a specified voliame, con- 
straints may be required. Thus, for example, in levitation 
experiments it may be extremely difficult to turn off the 
forces of constraint which keep the particles in a specified 
region. This means experiments which are sensitive to dis- 
turbances may be very difficult to perform if perturbation 
forces are required to be absent. A case in point are the 
drop dynamics experiments currently planned for Spacelab 
Mission 3. In these experiments an acoustic positioner is 
used to keep liquid drops in a field of view as well as for 
excitation of the droplets into oscillatory and rotational 
modes. However, after the droplets are excited via the 
acoustic field, it would be highly desirable in certain of the 
experiments to turn off the acoustic field in order to study 
the transition from one state of droplet motion to another. 
However, the acoustic field at the level of 10“4 g is 
required to keep the droplet in the field of view to cancel 
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the effects of spacecraft motion resulting from crew activity, 
drag deceleration , etc . 

In the analysis herein, it has been assumed that aero- 
dynamic drag on the particle is absent. However, aerodynamic 
drag will in most cases be present in the form of a Stokes- 
type of drag, given by 

D = 67rpvaV 
s 

where p is the density of the gas within which the particle 
resides, v is kinematic viscosity of this gas, a is the 
radius of the particle, and V is the velocity of the particle 
relative to the gas. If m denotes the mass of the particle 
and Pp the associated mass density, then we can rewrite the 
above as 



where t is a viscous time scale 

2 

_ 4 Pp^ 

^ 18 pv 

which provides a measure as to how long it takes for the 
particle to react to the effects of viscosity. Thus, if t 
is sufficiently large, viscous drag is unimportant, while if 
T is sufficiently small, viscous drag is important. Thus, we 
may conclude that viscous drag will be unimportant for cases 
in which particle mass density and radius are sufficiently 
large or in which gas density and kinematic viscosity are 
sufficiently small, with the reverse being true at the oppo- 
site extremes. 

The analysis herein has also neglected the effects of 
crew motion, venting of gases, and thruster firings. These 
effects can be included via specif ication in in the 
governing equations. 
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TIME (HRI 


FIG. I TYPICAL SEPARATION GROWTH BETWEEN A SPACECRAFT 
AND A PARTICLE IN LOW EARTH ORBIT. 
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A, INITIAL POINT 

B, END POINT 


FIGo I TYPICAL TRAJECTORY IH THE RELATIVE 
COORDINATE SYSTEM 
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DESIGN OF THE SPACE SHUTTLE DIGITAL AUTOPILOT 


AND RESULTING DYNAMIC ENVIRONMENT 

by 

Stanley Fay 

The Charles Stark Draper Laboratory, Inc. 


This compilation progesses through the development of the 
Space Shuttle phase plane controller from fundamental considerations, 
and it provides quantitative insight regarding the nature of the 
dynamic environment aboard the Shuttle. 

Figure 1 is a tutorial preliminary to the development of 
Euler's Equations. It introduces a vector torque as a time 
derivative of the angular momentum vector and then describes the 
necessary use of the Law of Coriolis to take the derivative. 

Figure 2 is a compact derivation of Euler's Equations. 

Figure 3 displays the assumptions used to simplify Euler's 
Equations for design of the Shuttle digital autopilot. 

Figure 4 displays further assumptions used in the DAP 
(digital autopilot) and presents a basic block diagram of the 
control loop. 

Figure 5 starts the basic development of the phase-plane 
controller of the DAP. It introduces the control variable U as 
torque per unit moment of inertia and demonstrates how phase- 
plane state drifts to the right for positive angular rate and 
drifts to the left for negative angular rate. 

Figure 6 derives the parabolic trajectory on the phase 
plane for the constant thrust of the attitude control jets. 

Figure 6 also demonstrates that one particular trajectory 
(dotted) is desired in order to achieve zero error in attitude 
position and attitude rate. 

Figure 7 introduces the concept of "switch lines" to dictate 
when to coast and when to reverse fire in order to ultimately 
approach zero error. 

Figure 8 introduces the concept of "deadband" and limit 
cycling, and how to slant the deadband switch lines in order 
to attenuate the limit cycle rate. 



Figure 9 shows how using parabolic deadband switch lines 
enable us to achieve the smallest limit cycle in one pass. 

Figure 10 and Figure 11 derive a common "one side of 
deadband" limit cycling which results from experiencing a 
constant torque disturbance. (Note that the larger the 
torque disturbance, the smaller the amplitude of limit cycle.) 

Figure 12 displays some numerical values of Shuttle 
limit cycle amplitude and period for various torque dis- 
turbances , T<j . 

Figure 13 presents an array of pertinent Shuttle facts. 

Figure 14 is a comparison of attitude stabilization 
accuracy of various systems. 

Figure 15 and Figure 16 are a derivation of the type and 
level of disturbance to be experienced aboard the Shuttle as 
a result of the firing of a single vernier jet for one 
minimum impulse. 

Figure 17 displays drag deceleration as a function of 
altitude . 

Figure 18 displays various accelerative g-levels due to 
various disturbance sources . 
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3-DIMENSIONAL SPACE 



FIGURE 1 



Mc=PiH = PbH + WibXH 
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FIGURE 2 
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SYSTEM DESIGN RATIONALE 


(1) IF THE CROSS PRODUCT TERMS, ARE SMALL COMPARED TO THE ACCELERATION TERMS, 

J3W3 >> {J2 — Ji) coi C02 

I.E., DIFFERENCES OF LARGE NUMBERS MULTIPLIED BY PRODUCT OF TWO SMALL 
NUMBERS, 

WE MAY WRITE: 

Tx ~ Jx 

Ty " '^y ‘^y DECOUPLED AXES 

Tz = Jz Wz 

(2) I F ANGULAR DISPLACEMENTS PER CALCULATION ARE KEPT SMALL, WE MAY SAY 
w T, AND THE SYSTEM CONTROL RELATES TORQUE TO ANGLE BY A SIMPLE 
DOUBLE INTEGRAL. NOTE THAT FOR LARGE ANGLES, AN ANGLE IS NOT A 
VECTOR. 

"NON-COMMUTATIVITY OF ANGLES." 


FIGURE 3 
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• FOR SMALL ANGLES, 0 = J //Td,2 

• ANGLES KEPT SMALL PER CALCULATION BY 12.5 Hz ITERATION RATE OF DAP 

• DRIFT OF ACTUAL ORIENTATION FROM DESIRED BECAUSE OF ASSUMPTIONS (1) 
AND (2) IS SENSED BY IMU AND NEW COMMAND ISSUED AT EACH ITERATION 


ACTUAL STATE 
— 


DESIRED STATE 
^ 


GUIDANCE 

ERROR 

RCS 

LAW 


ERRORS 


-►WE 


-►0E 


TO PHASE 
PLANE 


THRUST 



FIGURE 4 





PHASE PLANE 
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U = d = 03 


d = 03 


IF u = 0, NO TORQUE, 03 = 0 OR co = CONSTANT AND 0 INCREASES FOR POSITIVE 
oj AND DECREASES FOR NEGATIVE co 
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PHASE PLANE 



X 

J 


= U = CO 
0 = CO 


SLOPE 


dco 

d0 


ASSUME u- X 
J 


dco/dt _ 
d0/dt ~ 0 " w 

= CONSTANT FOR JETS: 


dco u 
d0 CO 


codco = ud0 

i (co2 - coj) = u (0 - 0 q) 



u = +1, COE must increase (HEAD UPWARDS) 
u = -1, COE MUST DECREASE (HEAD DOWNWARDS) 

COE IS +. 0E MUST INCREASE (HEAD TO RIGHT) 

COE 0E MUST DECREASE (HEAD TO LEFT) 

TO HIT 2ERO ERROR, DOTTED TRAJECTORY DESIRED. 

FIGURES 
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u = +1. WE INCREASES, MOVE UP PAGE 
u = -1, WE DECREASES, MOVE DOWN PAGE 
MOVEMENT ALWAYS CLOCKWISE 
NEVER HIT ZERO EXACTLY BECAUSE: 

1. ATTITUDE MEASUREMENT ERRORS 

2. RATE MEASUREMENT ERRORS 

3. MEASUREMENT DISTURBANCES OR NOISE 

4. JET QUANTIZATION FI RING TIME 


RESULT: 


1. CONTINUOUS LIMIT CYCLE AROUND ZERO 

2. FUEL WASTE 


FIGURE 7 
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ONE SIDE OF DEAD BAND LIMIT CYCLING 


JET IMPULSE: 


— TjAt = JAn (1) 

(JET FIRED IN AN ASSUMED NEGATIVE DIRECTION) 


DISTURBANCE TORQUE: +Tj = 


( 2 ) 


FROM (2): 


{DISTURBANCE TORQUE IS IN AN ASSUMED POSITIVE 
DIRECTION) 

dl2 Th 
dt J 


Th 

n-~ t + no 

j 

WHERE no * 1/2 An = - 1/2 


TjAt 


FROM (1) 


INTEGRATING (3): 


Jd t2 
J 

WHERE AT t = 0, 0 = 00 = ^DB 


0 = ^^2 *■ 



Dl FFERENTIATING (4) : 0 = t + 12o - 0 OR tmj„ - 


Td 2JTd 2Td 


TjAt 

PERIOD OF LIMIT CYCLE. T = -J— 

Td 


JET IMPULSE 
FIGURE 10 


(3) 


(4) 
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AMPLITUDE OF LIMIT CYCLE: 


FROM (4): A0 = = -1/2 t2 - ^Otmin 


Td T?At2 

= - 1/2 — + 


J 4T2 
d 


1/2 


= 1/8 


= 1/8 


T^^At2 

TdJ 

(JET IMPULSE)2 

TdJ 


NOTE: LARGER T^ YIELDS SMALLER Ad 
APPROXIMATE NUMERICAL VALUES: 

•=ivernler“25LB,At = 80MS 

MOMENT ARMS: ROLL, 10 FT 

PITCH, YAW, 40 FT 

ROLL Jx = 9 X 105 FT LB S2 
PITCH Jy = 72 X 105 FT LB S2 
YAW Jz = 75 X 105 FT LB S2 


FIGURE 11 



YAW AXIS 


GRAVITY GRADIENT 


Ad 

Td 

T 

1° 

0.006 FT LB 

3.7 H 

20 

0.003 FT LB 

7.4 H 

3.4 ^eS' 

6.5 FT LB 

12 S 



* TO LENGTHEN TIME BETWEEN JET FIRINGS, THE ACTUAL ADJUSTMENT IS FOR 
VA X DEAD BAND BY FIRING MORE THAN ONE JET. 


MIN0dB=O-1° 


FIGURE 12 
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SHUTTLE FACTS 

1. 38 MAIN JETS, 900 LB THRUST, TRANSLATION AND ATTITUDE 

^VERNIER JETS, 25 LB THRUST, ATTITUDE ONLY 

44 

la. APPROXIMATE MOMENT ARMS OF JETS FROM C.M. 

ROLL, 10 FT 
PITCH, YAW, 40 FT 

2. OMS ENGINES (TWO), GIMBALED, FIXED THRUST « 6,000 LB 

6° PITCH 
70 YAW 

3. ORBITER ANGULAR RATES: O.IO/STO 50 /S 

4. ATTITUDE DEAD BAND 0.lOTO5° 

5. ORBITER MASS APPROXIMATELY: 200,000 LB = g ggg sluqs 

9 

6. ORBITER MOMENTS OF INERTIA: 

ROLL (FORWARD). Jx 1.2 X 10® KG METER2 = 9 X 1q5 FT LB $2 
PITCH (RIGHT WING), Jy 9.7 X 10® = 72 X 105 

YAW (DOWN), J 2 10 X 10® = 75 X 10® 

FIGURE 13 
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ATTITUDE STABILIZATION 


1. SHUTTLE JETS: ±0.1° DEAD BAND 


2 . 

3. 


BOLT-DOWN CMC PACKAGE 


STABILIZED GIMBAL SYSTEM 


/ ACCURACY ± 

\ JITTER RATE ± 25 ®/S 

/ ACCURACY ± 1 
IJITTER RATE ± 1 ®/S 


EXAMPLE: APOLLO TELESCOPE MOUNT (ATM) ON SKYLAB 


- ATM TWO GIMBAL REQUIREMENTS: 


(VEHICLE AXIS) 
(PLATFORM AXES) 


- SKYLAB ATTITUDE VIA CMG 



POINTING 

(RAD) 

STABILITY 
(RAD/15 MIN) 

ROLL 

PITCH & YAW 

± 29 X 10-4 

± 12 X 10-4 
(250®) 

± 22 X 10-4 
± 26 X 10-4 



POINTING 

(RAD) 

STABILITY 
(RAD/15 MIN) 

ROLL 

PITCH & YAW 

± 29 X 10-4 

± 12 X 10-6 
(2.5®) 

±22X10-4 

±12X10-6 


FIGURE 14 
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ROLL AXIS: F = 25 LB 
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FIGURE 16 



DRAG DECELERATION 
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ACCELERATIVE g-LEVELS AND SOURCES 
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EXCEEDANCE STATISTICS OF ACCELERATIONS 


RESULTING FROM THRUSTER FIRINGS ON 
THE APOLLO-SOYUZ MISSION 
by 

George H . Fichtl 
Robert L. Holland 

NASA/George C. Marshall Space Flight Center 


INTRODUCTION 


Spacecraft acceleration resulting from firings of 
vernier control system thrusters is an important consider- 
ation in the design, planning, execution and post-flight 
analysis of laboratory experiments in space. In particular, 
scientists and technologists involved with the development of 
experiments to be performed in space in many instances 
required statistical information on the magnitude and rate of 
occurrence of spacecraft accelerations. Typically, these 
accelerations are stochastic in nature, so that it is useful 
to characterize these accelerations in statistical terms. 

This paper summarizes statistics of spacecraft accelerations 
associated with thruster firings that occurred during the 
Apollo-Soyuz mission. The motivation for performing this 
statistical analysis was simply to obtain an idea of the kind 
of exceedance statistics one might expect for vehicle 
accelerations associated with vernier control system thruster 
firings. To the best of our knowledge statistical data of 
spacecraft accelerations resulting from thruster firings are 
not available. The statistics of spacecraft accelerations 
for any given space mission depend on characteristics of the 
spacecraft (mass, control systems, etc.) mission profile, and 
a host of other parameters unique to the mission. 

Estimates of exceedance statistics of spacecraft 
accelerations can be obtained via Monte Carlo simulation of 
the mission, wherein the statistics of the forcing functions 
and other stochastic parameters germane to the mission enter 
the simulation via random selection from populations which 
possess statistical attributes of these forcing functions and 
parameters . 
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THRUSTER ACCELERATION TIME HISTORIES 


Acceleration at time t at a given point in a spacecraft 
as a result of thruster firings is given by 

g(t) = ^ X r + ^ X X r) + a (1) 


where is the rigid body angular velocity vector of the 
spacecraft, ? is position vector relative to the spacecraft 
center of mass, t is rectilinear acceleration, and (*) 
denotes time differentiation. During a thruster firing event 
^ and a are nonzero. After the event, these quantities vanish. 
Furthermore, the second term on the right side of equation 
(1) is usually two to three orders of magnitude smaller than 
the remaining terms (at least for the Apollo-Soyuz mission 
and typical Shuttle missions) . In addition, the typical 
thruster event duration time for attitude vernier control is 
less than a second. Thus, a time history of the magnitude of 
spacecraft acceleration, at a given point in the spacecraft, 
as a result of vernier control system thruster firings would 
be like a series of Dirac delta function-like acceleration 
spikes, corresponding to the thruster firings, with essen- 
tially zero acceleration between the spikes (see Figure 1) . 

Of course, the total time history would have additional con- 
tributions from crew activity, venting of gases and liquids, 
atmospheric drag, etc. This paper does not consider these 
additional contributions. 

Time histories of the kind depicted in Figure 1 were 
calculated for specific periods of time during the Apollo- 
Soyuz Mission for postflight analysis of the experiments per- 
formed on that mission. The calculations involved the 
estimation of w and a via the calculation of the net force 
and net torque associated with each thruster firing event. 

This was possible because the location and reaction force of 
each thruster relative to the center of mass was known. A 
discussion of the Apollo-Soyuz vernier control system and 
thrusters is provided in the reference cited at the end of this paper 

Four acceleration time histories were selected for 
analysis, namely those associated with experiments MA-044, 

085, 041 and the science demonstration (Sci-Dem) . These runs 
were approximately 3 to 17 hours in duration. The thruster 
firing time for each thruster firing event was 0.015 or 
0.031 sec. 
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STATISTICAL CALCULATIONS 


The quantity g is a vector. We have performed statisti- 
cal calculations on the individual components and the magni- 
tude of §. In this section we present our statistical 
calculations on the magnitude of We seek the answer to 
two questions. First/ what is the average rate of exceedance 
of vehicle acceleration at some assigned level g^? (See 
Figure 1) . Second, what is the statistical distribution of 
the time interval between adjacent acceleration events? 

The answer to the first question would provide an 
experimenter an estimate of the average number of times a 
critical acceleration level would be exceeded during an 
experiment with duration time t@. Furthermore, by making 
assxmptions about the distribution of the number of times a 
critical acceleration level is exceeded in time tg an esti- 
mate can be made of the risk of exceeding a critical acceler- 
ation level during the experiment. 

The answer to the second question would provide the 
experimenter with an estimate of the probability of the time 
interval between acceleration events being less than or equal 
to some critical time. 


ACCELERATION EXCEEDANCE STATISTICS 


We define the exceedance rate N of the magnitude of the 
acceleration vector, g say, to be the average number of 
acceleration events that exceed level gj^ per unit time. To 
calculate this quantity from the data we merely count the 
number of events that exceed level gj^. In an attempt to 
combine the exceedance calculation for the four time histo- 
ries that were examined we (1) scaled the exceedance rate N 
with the standard deviation, of the time interval between 

acceleration events (tx, t2, ta, ... in Figure 1) and (2) sub 
tracted the mean value of g, g say, from g and scaled the 
resulting difference with the standard deviation of g, o„ say 
Thus, Figure 2 provides plots of results of calculations^of 
nondimens ional exceedance rate Na.p as a faction of the non- 
dimensional unbiased accelerations, (g - g)/Ug. 
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DISTRIBUTION OF TIME INTERVALS 


To calculate the distribution of time interval between 
acceleration events, T say, straightforward counting pro- 
cedures were used. Here again in an attempt to combine the 
results of_our calculations we subtracted the mean time 
interval, T say, from the time interval and scaled the 
resulting difference with Ciji. Thus, Figure 3 provides plots 
of the distribution function of unbiased time interval 
between acceleration events, (T - T)/aiji. 


DISCUSSION 


It appears that some order can be brought about in the 
statistics of spacecraft vernier control thruster firing 
statistics via the scaling introduced in the previous section. 
This is especially true for the thruster time interval dis- 
tribution function. The thruster acceleration exceedances do 
not collapse as weljL as the time interval distribution 
functions for (g - g)/ag < - 1. However, this may not be 
important because the important part of the curve is for 
(9 “ 9)/<^g > 0,_i.e., the part for g > g. The part associ- 
ated with (g - g)/ag < 0 encompasses 0 £ g < g. 


We do not claim these results to be universal because of 
variation in mass, crew activity, control system parameters, 
etc., from one spacecraft to the next. If results like these 
could be derived for planned missions, they could play a 
useful role in the planning of experiments which are sensi- 
tive to spacecraft dynamics. These statistics permit the 
calculation of the risk associated with g > g^ for a given 
experiment duration time tg, where g^ is a critical value of 
g to the experiment. To do this one needs to assume that the 
exceedances of the g process above level g^ arrive inde- 
pendently. We note by Q(tg) the number of exceedances of g 
at level g^ over the experiment duration time tg. Clearly, 
the process Q(tg) is a Poisson process, and the probability 
of Q(tg) being less than or equal to an assigned value (for 
example q) is given by 


(Xt -At 

P(Q(tg) < q , tg) = e ® 


( 2 ) 


where X is a parameter. The probability of no exceedance of 
the g process above the critical value g^ in time interval tg 
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follows by setting g = 0 in the above equation, so that 


P(Q) (tg) = 0 , tg) = e ® (3) 

By definition of the Poisson process we set 
X = N 

Now the risk R that the g process will exceed the critical 
value of g^ at least once during an experiment of duration 
time tg is thus 


-Nt 

R = 1 - e ® (4) 

Upon specifications of risk R and tg the value of N can be 
calculated. A curve like that shown in Figure 2 can be used 
to calculate g upon specification of Ciji, g, and Cg. These 

latter parameters as well as the Na^ versus (g - ^)/<7g curve 
will in all likelihood vary from mission to mission of a 
given spacecraft and possibly vary within a mission. A curve 
like that shown in Figure 3 can be used to estimate the 
distribution of time intervals between g-events upon speci- 
fication of ^ and 0 ^. This distribution would provide an 
estimate of extreme values and typical values of T. Thus, 
the scientist obtains a statistical description of the 
dynamic environment in terms of the risk of exceeding a 
critical g-level and the distribution of time intervals 
between g-events. 


REFERENCE 


Apollo Test Project Operation Data Book, Vol. 
Lyndon B. Johnson Space Center, Houston, 


1, MSC 07765, 
Texas, 1973. 
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FIGURE 1 CALCULATION OF THRUSTER FIRING STATISTICS 



(oas) (6g0l) 



FIGURE 2 ASTP THRUSTER ACCELERATION EXCEEDANCES 













FIGURE 3 ASTP THRUSTER TIME INTERVAL DISTRIBUTION FUNCTION 




ELASTIC BODY DYNAMICS 


by 

Ben W. Holder 

NASA/Lyndon B. Johnson Space Flight Center 


Most of the structural dynamics resources allocated to 
the Space Shuttle have been concentrated on the flight events 
which result in critical structural loads and/or minimum 
control stability margins. Since these events are primarily 
sxab-orbital , the data base of interest to those involved in 
orbital experimentation is somewhat limited. A brief dis- 
cussion of available data is given in the following para- 
graphs. 

Three types of on-orbit thrusting are planned. These 
involve the reaction control system (RCS ) , vernier reaction 
control system (VRCS) , and the orbital maneuvering system 
(QMS) . Of the three, the OMS provides the largest acceler- 
ation levels and in all likelihood, the most severe elastic 
responses. However, the OMS is used primarily for orbital 
insertion and de-orbit maneuvers and, therefore, is not likely 
to be a major source of disturbance during on-orbit experi- 
mentation. The RCS was designed as the primary system to be 
used for on-orbit control. However, present planning calls 
for use of the VRCS rather extensively for this purpose since 
it is more efficient from a fuel usage standpoint. The VRCS 
also provides small control system deadbands and lower 
acceleration levels. A disadvantage of this system is its 
lack of redundancy. Loss of a thruster can render the system 
ineffective. For this reason commitments to use only the 
VRCS for protracted periods has been discouraged. Other on- 
orbit disturbance sources which can produce elastic responses 
are on-board mass movements, most of which are associated 
with crew motion or remote manipulator system activity. 

Two major structural configurations will exist while on- 
orbit: Orbiter payload bay doors open and Orbiter payload 

bay doors closed. Up-to-date structural math models are 
easily developed for the latter case since it represents 
minor changes to models that are required for response 
analyses for several mission phases. An up-to-date open door 
structural dynamic math model does not presently exist. How- 
ever, a model which was developed for studies concerned with 
the Spacelab IPS is available for preliminary analysis 
purposes . 
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The development of the on-orbit elastic response environ- 
ment requires a detailed knowledge of the disturbance inputs 
as well as the structural math models and disturbance sources, 
some of which were discussed earlier. As far as the author 
can determine, data describing crew motions and thrusting 
duty cycles for the Orbiter are not presently available. 
Certain basic information can be derived from available data 
sources which should be of interest to those involved in on- 
orbit experimentation. First, the rigid body acceleration 
levels due to QMS, RCS, and VRCS firings are given in Figures 
1 and 2. These levels contain an allowance for elastic 
response. Second, Figures 3 and 4 provide examples of trans- 
fer function information for the closed-door math model. The 
transfer functions show the ratio of the acceleration response 
in the Orbiter mid-bay to steady-state sinusoidal force input 
at the various thruster locations. These data are useful in 
determining the frequency content of the response. Similar 
data exist for the open-door math model. 

In summary, information regarding structural responses 
during on-orbit operations is limited. Although estimates of 
peak acceleration levels and the associated frequency spectrum 
in the payload bay due to thrusting of the various control 
system thrusters have been made, the actual levels and time 
histories must be based on updated structural math models and 
a detailed knowledge of the input forcing functions. 
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SIGN! FICANT ON-ORBIT FORCING FUNCTIONS 


• ORBITAL MANEUVERING SYSTEM 

LOAD FACTOR i = EXTERNAL FORCE IN i DIRECTION 

VEHICLE WEIGHT 


CARGO MAXIMUM LIMIT LOAD FACTORS/ANGULAR 
ACCELERATIONS FOR OMS OPERATION 



LOAD FACTOR 

ANGULAR ACCEL. 
RAD/SEc2 

CARGO WEIGHT 

Nx 

Ny 

Nz 

•• 

Ox 


h 

<32K 

-0.27 

±0.005 

-0.09 

±0.01 

±0.005 

±0.005 

>32K 

-0.24 

±0.004 

-0.07 

±0.01 

±0.005 



FIGURE 1 






















• REACTION CONTROL SYSTEM 


CARGO MAXIMUM LIMIT LOAD FACTORS/ANGULAR 
ACCELERATIONS FOR RCSA/RCS OPERATION 


SYSTEM 

LOAD FACTOR 

ANGULAR ACCEL. 
RAD/SEC2 

Nx 

Ny 

Nz 

•• 

&v 


RCS (870 LBS) 

0.019 

±0.022 

m 

±0.021 

M 

±0.014 


-0.016 


HB 




VRCS (25 LBS) 

+0 


+0.0002 

±0.0007 


±0.0003 


1 

o 

■1 

-0.0 





CENTER-0 F-M ASS AT (1080, 0, 400) 

ANGULAR ACCELERATIONS MAY OCCUR SIMULTANEOUSLY 


FIGURE 2 
























EMPTY M54D ATTACH POINT ACCEL TF 
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FREQUENCY (HZ) 


X RESPONSE AT X = 1002.13.y = -94.. z = 414. 

TWO SINUSOIDAL THRUST FORCES WITH UNIT 
AMPLITUDE AT OMS POD CG ALONG X-AXIS 


FIGURE 3 
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EMPTY M54D ATTACH POINT ACCEL TF 



FREQUENCY (HZ) 


Z RESPONSE AT x = 1002. 13. y = 94. .z = 414. 

TWO SINUSOIDAL FORCES WITH UNIT AMPLITUDE 
AT TWO NOSE GEAR TRUNNIONS ALONG Z-AXIS 


FIGURE 4 
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CREW ACTIVITIES IN SPACE 


by 

Lt. Col. Guian S. Bluford, Jr. 
NASA/Lyndon B. Johnson Space Flight Center 


INTRODUCTION 


One of the mission requirements of the Space Shuttle is 
to serve as a working platform for experiments in space. 

Many of these experiments will be performed by crewmembers 
(mission specialists and payload specialists) in a general- 
purpose laboratory called Spacelab. Part of Spacelab con- 
sists of a pressurized module, accessible from the Orbiter 
cabin through a transfer tunnel, which will provide a shirt- 
sleeve environment for performing experiments in space. All 
nonexperiment-related activities or housekeeping activities will 
be done in the Orbiter, while most of the mission-related 
activities (experiments) will be done in Spacelab. 

In order for experimenters to design their experiments 
to best utilize the capabilities of the Orbiter, the Spacelab, 
and the crew, it is necessary that they understand the work- 
ing environment in the Orbiter and in Spacelab. In addition, 
the experimenters must be knowledgeable of the housekeeping 
activities required of the crew. 


ORBITER 


The Orbiter is the crew and payload carrying section of 
the Space Shuttle. It is a delta-winged spacecraft/aircraft, 
about the size of a DC-9 airliner. The forward part of the 
Orbiter contains the pressurized section in which the crew 
works, while the aft five-sixths of the Orbiter consists of a 
payload bay in which the Spacelab module will sit (see 
Figure 1) . 

During launch and entry, the acceleration stresses on 
the crew and payloads will never exceed three g's. On orbit 
the acceleration levels will be much lovjer ( <.06 g ) and 
will be based on mission requirements and Orbiter capability. 

The Orbiter will also be subject to random vibrations on 
its exterior surfaces, particularly during launch. These 
vibrations will be generated by acoustic noise (generated by 
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the engine exhaust) and by aerodynamic noise (generated by 
airflow) . These fluctuating pressure loads will cause 
structural vibrations which will be transmitted to the mid- 
fuselage payload section of the Orbiter and eventually to the 
payload. The magnitude of these vibrations will be based on 
the transmission characteristics of the Orbiter payload 
support structure and the interactions between each payload's 
mass, stiffness, and center of gravity. The estimated 
values of thermal, pressure, contamination, and electro- 
magnetic environments in the payload bay for various aspects 
of flight are given in Space Shuttle System Payload Accommo- 
dations, JSC-07700, Voliime XIV, for user reference. 

The pressurized portion of the Orbiter provides separate 
working and living quarters for the crew with ordinary air 
(21% oxygen, 79% nitrogen) provided at standard sea-level 
pressures of 14.7 psi and temperatures of 11°G to 27°C 
(Figure 2). The hiamidity is controlled, and odors and carbon 
dioxide are continuously filtered out. 

The upper section of the cabin is the flight deck, fromwhich 
the Shuttle is controlled and most external payloads are handled. 
During on-orbit operations, the crew monitors and controls 
Shuttle systems, performs rendezvous and stationkeeping, 
operates the payload bay doors and radiator systems, and 
manages the Orbiter 's interface with payloads and Spacelab. 

Below the flight deck is the crew living area or mid-deck. 

In this mid-deck area are the sleeping stations, the galley, 
and the washroom. From the back of the mid-deck is a hatch 
which leads to the tunnel for entry into the pressurized 
module of Spacelab. 


SPACELAB 


Spacelab consists of a module and pallet sections used 
in various configurations to suit the needs of a particular 
mission (Figure 3) . The pressurized module provides a 
shirtsleeve working environment for four crewmembers. The 
Spacelab subsystem equipment, such as environmental control 
system, power systems, data-handling systems, etc., are 
located in the forward end of the module, leaving about 60% 
of the volume available for experiments. In addition, some 
mission— dependent equipment, such as a top airlock and an 
optical window/viewport assembly can also be flown. More 
detailed information on accommodations can be found in the 
Spacelab Payload Accommodations Handbook (European Space 
Agency SLP/2104) . 
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The pallets accoramodate experiment equipment for direct 
exposure to space. Each standard pallet segment is 3 meters 
long, and two or three pallets can be connected together to 
form a single pallet train. When no module is used with the 
pallets, a cylindrical "igloo” mounted on the end of the for- 
ward pallet provides a controlled, pressurized environment 
for Spacelab subsystems normally earried in the module. A 
mission-dependent instrument pointing subsystem can be 
attached to a pallet for precision payload instrument point- 
ing. Pallet equipment, without a Spacelab module, is 
operated remotely from the Orbiter aft flight deck or from 
the ground. 


CREW OVERHEAD ACTIVITIES 


Besides performing mission-related experiments in space, 
the crew will have collateral duties to be performed in the 
Orbiter. These activities include sleeping, eating, exer- 
cising, commionicating with the ground, and Orbiter systems 
management. These activities must be integrated into the 
daily routine or timeline of the crew. A synopsis of these 
overhead activity requirements is given as follows: 


Pre-sleep/Post-sleep Activities 


Approximately 45 minutes will be allotted for each pre- 
sleep and post-sleep activity period. These activities 
include stowing or unstowing equipment, reviewing checklists 
and crew activity plans, and personal hygiene. 


Sleep 


Eight hours per day of uninterrupted sleep will be 
allotted to each crewman. The sleep stations in the mid-deck 
of the Orbiter can accommodate four crewmembers simultane- 
ously. 


Meals 

One hour per meal (three meals per day) will be allotted 
to each crewmember for food preparation and dining. Four 
crewmembers can be accommodated simultaneously in the galley 
area. 


Exercise 

A limited amount of physical exercise is 
anticipated for each crewman on flights exceeding eight 


183 


j 



days. Exercise will be performed on a portable treadmill or 
isokinetic exerciser. 

Communications/Handover (Two-Shift Operations Only) 

Approximately 15 minutes per day will be allotted to 
reporting Orbiter systems and cons\amables status to the 
ground. The crew required for this status report will be the 
commander and/or pilot and a mission specialist. 

For two-shift operations, 15 minutes per day will be 
required for hand-over operations on which the old shift will 
brief the new shift on spacecraft status, continuing activi- 
ties, and any changes in mission operations to be performed 
during the new shift. 

Orbiter Systems Management 

Most Orbiter systems management functions will be per- 
formed by the commander and/or pilot with some collateral duties 
being performed by a mission specialist . These collateral duties 
include fuel cell purging (5 minutes in the auto mode or 25 
minutes in the manual mode every 8 hours) , CO 2 absorber 
changeout (5 minutes every 12 hours for four crewmen) , water 
dumps (as scheduled) , inertial measuring unit alignments 
(15 minutes every 8 hours) , trash management (as required) , 
and housekeeping duties (at least one hour per day) . 


MISSION PLANNING 


In preparing a crew activity plan for payload activities 
on a Spacelab mission, it should be anticipated that crew 
work efficiency will initially be low but will rapidly 
increase as the crew adjusts to the zero-g environment. Pay- 
load specialists will be dedicated full time on Spacelab 
experiments, while mission specialists will be dedicated 
approximately 80% of the time. Orbiter activities and 
detailed test objectives will be integrated with payload 
activities. 
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FIGURE 2 
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ACOUSTICS AND VIBRATION ENVIRONMENTS FOR 
LABORATORY EXPERIMENTS IN SPACE 
by 

W. P. Rader 

Martin Marietta Corporation 
Denver, Colorado 


INTRODUCTION 


This paper briefly desGribes some of the basic param- 
eters of random noise and vibration, and presents typical 
environments for the launch phase and orbital operations. 

For the latter, both acoustically induced and structure-borne, 
thruster-induced vibration are addressed, using data obtained 
during the Skylab and Titan programs. 


LAUNCH ENVIRONMENTS 


In order for experiments to perform the scientific 
functions for which they were designed, the equipment must 
first survive the launch/boost environments. Typical 
acoustic and random vibration criteria [1] for Spacelab 
equipment are presented in Figure 1. In Figure la, the 
acoustic levels internal to the payload bay are applicable to 
pallet-mounted experiments. The noise reduction of the 
Spacelab module attenuates the environment to the lower level 
shown {Figure la, curve B) . Typical predicted random vibra- 
tion criteria for pallet- and module-niounted equipment are 
shown in Figure lb. Note that the pallet specification is 
applicable for equipment mounted at pallet hard points. For 
experiments with large, lightweight surface areas, the 
response of the experiment to direct acoustic excitation must 
be considered also, since the vibration levels produced may 
be higher than those transmitted from the pallet. 


ORBITAL ENVIRONMENTS 


The primary emphasis in this paper is on the acoustic 
and vibration environments which may interfere with experi- 
ment function or scientific data acquisition during opera- 
tions in space. Very limited measured data are available. 
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since, in the past, the emphasis for conventional launch 
vehicles has been on defining the severe launch environments, 
and the data acquisition systems used had inadequate dynamic 
range to measure the low-level orbital environments. The 
reusable Shuttle, with suitable instrvimentation systems and 
on-board recording capability, provides the potential to 
alleviate the problems associated with the limitation of 
telemetry systems. 

Data obtained in the Multiple Docking Adapter (MDA) 
during the Skylab Missions have been selected as typical of 
the acoustical environment which may be expected in Spacelab. 
As points of reference, the overall sound pressure levels in 
the MDA ranged from roughly equivalent to those in a noisy 
office or department store to those on a noisy, busy street. 
The noise-producing equipment in Skylab included fans and 
pumps associated with the refrigeration and environmental 
control systems, the rate gyro package, and gears and stepper 
motors associated with experiments. 

Unfortunately, no random vibration measurements were 
obtained during Skylab orbital operations. Therefore, esti- 
mates of the vibration levels have been made based on ground 
test data and adjustments for differences in acoustic levels 
measured during the ground test and during the mission. The 
general technique is often used in predicting random vibra- 
tion environments for components on new programs. The appli- 
cable data were obtained from References 2 and 3, and the 
results are shown in Figure 2 . 


THRUSTER- INDUCED VIBRATION 


The response of spacecraft to thruster actuation is 
potentially more detrimental to science experiments than is 
the equipment noise. Propulsion systems used for orbital/ 
attitude control typically range in thrust capability from a 
few pounds to several thousand pounds. One might expect the 
effect to be limited to a steady-state acceleration; however, 
the actuation of valves and of pressure fluctuations in the 
propulsion process produces a dynamic forcing function. The 
effect is illustrated by an example using data from a Titan 
transtage . 

Data were obtained from several low-frequency (0 to 
50 Hz) accelerometers mounted on primary airframe structure 
during the transtage engine firing. A typical vibration time 
history is shown in Figure 3. The resulting vibration 
spectrum is presented in Figure 4. Since experimenters are 
often concerned with displacements rather than accelerations. 
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the maximxam displacement of the vibration has been estimated. 
Note that the vibration spectrum has a significant peak at 
approximately 8 Hz. Referring to the acceleration time 
history (Figure 3) , the maximum acceleration amplitude at 
this frequency is approximately .05 g peak and the associated 
maximum displacement is approximately 0.0008 in., 0 to peak. 


CONCLUSIONS 


Some of the basic parameters associated with fandom 
noise and vibration have been briefly discussed, and typical 
acoustic and vibration environments for both the launch phase 
and orbital operations have been presented. Too little infor- 
mation is available on the orbital environments. Clearly, 
measurement systems to acquire these data are needed on early 
shuttle flights, as well as improved data analysis techniques 
to determine sources of excitation and transmission of noise 
as a function of frequency. 

It is hoped that the information presented will provide 
experimenters with a clearer understanding of acoustics and 
random vibration. Clearly, we dynamicists need to develop a 
better understanding of the experiments and interaction with 
the environments in order to recognize what may be detri- 
mental and how to alleviate the problems. 
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191 


4 






FIGURE 1b. launch VIBRATION LEVELS 


FIGURE 1. SPACELAB LAUNCH ACOUSTIC AND VIBRATION ENVIRONMENTS. 
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FIGURE 4 VIBRATION RESPONSE TO THRUSTER FIRING 


196 


i 



APPENDIX A 

AGENDA FOR WORKSHOP ON SPACECRAFT DYNAMICS 
AS RELATED TO LABORATORY EXPERIMENTS IN SPACE 

May 1, 1979 

8:00 Registration 

8:30 Welcome .......... W. W. Vaughan 

8:35 Workshop Objectives 

and Plan . G. H. Fichtl 

MORNING SESSION Flight Experiments 
Chairperson: B. Antar 

8:55-9:00 Introduction B. Antar 

9:00-9:15 Critical Point Experiment . M. Moldover 

9:15-9:30 Liquid Helium 

Experiments P. Mason 

9:30-9:45 Geophysical Fluid 

Dynamics W. Fowl is 

9:45-10:00 Drop Dynamics ....... D. Elleman 

10:00-10:15 Atmospheric Cloud 

Physics J. Anderson 

10:15-10:30 Break 

10:30-10:45 Combustion Processes ... A. Berlad 

10:45-11:00 Boiling Heat Transfer . . . W. Frost 

11:00-11:15 Tribology Experiment- ... R. Gause 

11:15-11:30 Materials Processing 

in Space ........ R. Naumann 

11:30-12:15 Convection in Fluids at 

Reduced Gravity ..... S. Ostrach 

12:15-1:15 Lunch 
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AFTERNOON SESSION Spacecraft Dynamics and Low Gravity 

Fluid Mechanics 

Chairperson: F. Collins 

1:15-1:20 Introduction . F. Collins 

1:20-1:40 Rotational Fluid 

Mechanics R. Cans 

1:40-2:00 Vehicle Motion Induced-g 

Effects on Fluids .... R. Dressier 

2:00-2:10 Orbital Density 

Environment R. Smith 

2; 10-2: 40 Orbital Mechanics L. Mullins 

2:40-3:00 Relative Motions Inside 

Spacecraft B. S. Perrine 

3:00-3:15 Break 

3:15-4:00 Spacecraft Rigid Body 

Dynamics S. Fay 

4:00-4:30 Elastic Body Dynamics ... B. W. Holder 

4:30-4:50 Crew Activity F. Bluford 

4:50-5:05 Measurement of Spacecraft 

Dynamics S. Fay 

5:05-5:30 Acoustics and Vibration . . P. Radar 

5:30 Adjourn 

6:00 Social . . Redstone Arsenal Officers Club 

May 2, 1979 

MORNING SESSION 

Chairperson: G. H. Fiehtl 

9:00-9:05 Introduction G. H. Fiehtl 

9 : 05-9 : 3 0 Miniature Electrostatic 

Accelerometer W. Lang 

9:30—9:50 Gravity Wave Detection . . R. Drever 
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9:50-10:00 


E. Eller 


ORI Activities in Materials 
Processing in Space . . . 

10:00-10:30 Summary of Previous 

Day's Meeting G 

10:30-12:30 Group Discussion 

12:30 Adjourn 
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